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Abstract 

The  conmunication  capacity  of  a single  wideband  satellite 
channel  can  be  simultaneously  shared  by  a number  of  users  by  means  of 
spread-spectrum  multiple-access  (SSMA) . In  phase-coded  SSMA  the  multiple- 
access  capability  is  provided  by  phase  modulating  a distinct  signature 
sequence  onto  the  user1 3 carrier  which  spreads  the  user  data  over  a 
wide  bandwidth.  All  of  the  Important  code  parameters  for  the  analysis 
of  such  a system  can  be  derived  from  the  aperiodic  correlation  functions 
of  the  signature  sequences.  The  asymptotic  behavior  of  such  code  param- 
eters is  considered  for  random  binary  sequences  for  which  the  sequence 
length  grows  very  large.  New  sets  of  pseudo-random  or  m-sequeuces  with 
optimal  aperiodic  autocorrelation  and  cross -correlation  properties  are 
obtained.  The  relationship  between  the  first  few  central  moments  of  the 
aperiodic  correlation  functions  and  the  characteristic  polynomials 
generating  the  m-sequences  is  analyzed  and  the  results  are  compared  with 
actual  sequence  data. 

Gauss'  product  of  cyclctomic  co  . u5  ujed  to  establish  new 
analytical  results  on  the  periodic  correlation  properties  of  Gold  sequences 
and  Kasaml  sequences  yielding  subsets  of  sequences  whose  correlation 
parameters  satisfy  tighter  bounds  than  previously  established  for  the 
entire  sequence  sets.  Numerical  data  on  the  relevant  correlation  parameters 
is  obtained  for  a large  number  of  good  signature  sequences. 
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CHAPTER  1 
INTRODUCTION 

1.1.  Spread-3pectrum  multiple-access  cocmunication 

In  recent  years  there  has  been  a large  Increase  in  the  number 
of  satellite  conrounication  programs  for  civilian  as  well  as  military 
purposes.  Pritchard  (1977)  lists  32  satellite  systems  currently  in  use. 

A major  advantage  of  a synchronous  or  geostationary  satellite  is  the 
wide  coverage  area  which  enables  a number  of  widely  dispersed  users  to 
have  simultaneous  access  to  the  satellite  transporter.  The  user 
channel  separation  can  be  achieved  in  a variety  of  ways: 

1)  Frequency-division  multiple-access  (FDMA)  is  a common  form  of  multiple- 
access  where  each  ground  station  has  a different,  precisely  determined 
carrier  frequency.  Siugle-channel-per-carrier  as  well  as  multiplexing 

to  multi-ch&nne Is -per -carrier  is  possible  and  fits  wall  within  non-digital 
terrestial  communication  networks. 

2)  If  precise  time  cooperation  between  the  transmitting  stations  is  possible 
one  can  adopt  the  very  efficient  time-division  multiple -access  (TDJiA)  in 
which  each  user  has  the  same  carrier  but  operates  in  a different  time 

a lot.  TDKA  fits  better  in  digitised  coccrunication  networks. 

3)  In  Code-division  multi pie -access  the  channel  separation  is  primarily 
due  to  coding  while  no  precise  frequency  or  timing  cooperation  between 
the  transmitting  stations  is  necessary.  Applications  are,  among  others, 
in  tracking  and  data-relay  systems  (Staapfl,  et  al.,  1970),  air  traffic 
control  (Stiglits,  1973)  and  military  satellite  cocouni cation  systems 
(Gerhardt,  1973).  Spread-spcctrus  techniques  are  characterised  by  the  use 
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of  a great  deal  more  modulated  RF  bandwidth  than  normally  would  be  required 
for  transmitting  the  user  Information.  The  two  most  common  forms  of 
spread-spectrum  techniques  employed  in  CDMA  are  frequency-hopping  and  direct 
sequence  modulation  (Dixon,  1976).  Frequency -hopping  was  used  in  the  TATS 
modulation  system  and  is  described  in  Drouilhet  (1969).  In  direct 
sequence  modulation  the  multiple-access  capability  is  provided  by  a high- 
rate  code  sequence  which  phase  modulates  — together  with  the  data 
sequence  — the  user's  carrier.  The  receiver  station  will  recover  the 
information  by  means  of  correlation  techniques.  The  phase-coded  spread- 
spectrum  multiple-access  (SSMA)  method  is  very  attractive  for  communication 
systems  which  also  require  protection  against  malicious  interference  and 
unauthorized  listening  (Gerhardt,  1973). 

The  performance  of  the  SSMA  system  depends  on  the  correlation 
properties  of  the  high-speed  signature  sequences.  In  the  past  (Aein,  1965), 
(Blasbalg,  1965)  most  of  the  attention  was  focused  on  the  periodic 
correlation  properties  of  the  signature  sequences  because  in  most  cases 
a synchronous  communication  model  was  assumed.  In  the  next  section  ve 
Introduce  an  asynchronous  SSMA  system  model  earlier  presented  by 
Purs  ley  (1974)  --  which  makes  it  possible  to  identify  the  correlation  param- 
eters of  interest  for  the  communication  performance  as  well  as  for  the 
synchronization  performance. 

1.2.  Ehase -coded  SSMA  system  model 

Ve  will  consider  the  SSMA  system  model  as  given  in  Figure  1 for 
K transmitting  stations  or  users.  The  i-th  user's  data  signal  b,(t)  Is  a 
sequence  of  positive  and  negative  pulses  of  duration  T and  unit  amplitude. 


pcctrura  multiple-access  system  model. 


bi(t)  = JL  Vi  px<“«) 


(1.1) 


where  b^  ^ €{+l,  -1}  denotes  the  i-th  user's  information  bit  stream  and 
pulse  pT(t)  * 1 for  0 — t < t and  PT(fc)  a 0,  otherwise.  Hie  code  wave 
form  ^iich  binary  phase  modulates  the  user's  RF  carrier  frequency  can  be 
expressed  as 


a (t)  - E u (i)p  (t-jTe) 

jsa-ao  J 1c 


(1.2) 


where  {u,^}  represents  the  discrete  signature  sequence  of  the  i-th  user, 
and  has  period  p * T/Tc  and  elements  of  {+1,-1}. 

For  an  asynchronous  system  where  no  timing  reference  for  the  K 
users  is  assumed,  the  received  signal  at  the  message  destination  can  be 
expressed  >is 


K £ 

r(t)  = E (2P)  a (t-T  )b  (t-T  )co8(m  t +<t>  ) + n(t)  . 

1 3.  1 1 Cl 


(1.3) 


Here  n(t)  represents  the  channel  noise  which  we  assume  to  be  a white 
Gaussian  process  with  two-sided  spectral  density  NQ/2,  iuc  represents 
the  common  RF  canter  frequency  and  P the  common  signal  power.  Unequal 
signal  powers  can  easily  be  incorporated  in  the  results.  If  the  received 
signal  r(t)  is  the  input  to  a synchronised  correlation  receiver  matched 
to  the  i-th  user  signal,  the  output  at  sample  moment  t * T is  given  by 
Purs  ley  (1974)  as 

I K T 

Zi^  = (l)  ^bi,0  * T + 2 Si,k^Tk^  + J*  n(t)a  (t)cos  ^t  dt  . 

k“l  * 0 


(1.4) 


where  Si  ^(t^)  — with  [u^3  ■ {u^1^}  and  £ v ^ } - {l^^}  — equal, 
f°rtkIcSTkS<1k  + l)Ie 

'\,Ol>u>v<VIc  + “t/wl1  ' 'v¥'TrtV>“  % 


Si,k(V  * < 


lf  bk,0  * \,-l 

bk,Ot5u,v(VTc  + [®u,»<1k+l)  * V1(Vl(V‘k:c)lc“  ,l 


L 


1£  »k,0  * \,-l  • (l-5) 


The  discrete  cross -correlation  functions  0 (X)  and  8 (X)  for  the 

u,v%  u,vN 

sequences  u a {u^}  and  v = {v^}  are  defined  as 


.d  (i)  =*  C (X-p)  + C (X)  , 0-^  X ^ p-1 

J u,v'  u,v'  r/  u,vv  * r 

'§  (X)  «=>  C (X-p)  - C (X)  , 0 ^ X 5 p-l 

u,vv  7 u,vv  u,vs  7 * r 


(1.6) 


where  C (X)  denotes  the  aperiodic  cross -correlation  function  defined  as 
u,v'  7 r 


r 


C (X)  « \ 

u,vv  ’ 


p-l-X 

A "jv* 


p-14X 

z Vxvi 
j-o  j * j 


OS  ^ p-l 

1-p  5 i < 0 

ix|  * p 


(1.7) 


Observe  that  y(X)  » Cv  u(-X),  thus  the  periodic  (or  even) 

cross -correlation  function  9^  v(X)  satisfies  9^  V(X)  “ 8y  u(p*X)  'diere 

as  function  9 (X)  satisfies  9 (X)  =»  - 0 (p-X)  --  hence,  Massey  and 

u,v  u,v  v,u'r  7 

UUran  (1974)  called  the  latter  the  odd  cross -correlation  function. 
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In  order  for  the  correlation  receiver  to  operate  properly  it 
should  be  synchronized  with  the  i-th  user's  signal.  The  synchronization 
decisions  are  derived  from  Z^(T)  in  (1.4)  which  will  contain  the  periodic 
autocorrelation  function  0 (2)  * 0 (2)  or  the  odd  autocorrelation  function 

U VI  }U 

0^(2)  * 0u  u(2)  as  long  as  the  receiver  is  not  yet  synchronized  or  has  lost 
synchronization.  The  autocorrelation  functions  will  also  appear  in  Z^(T) 
when  multipath  interference  is  present  in  the  channel  (Massey  and  Ohran,  1969). 

1.3.  Worst-case  performance 

The  function  S.  , (t  ) in  equation  (1.4)  achieves  a maximum  value 
with  respect  tx>  and  q>k  whenever  » i^T,, , for  some  integer  and  q>k  “ 0. 
The  resulting  values  of  ^(t^)  are  + P^  °r  - v^k^*  Kence  t^Je 

maximum  value  ot  the  error  probability  Pr{Zi(T)  > 0 / b^  Q = - 1}  will  be 

minimized  by  selecting  a set  of  signature  sequences  for  which  the  peak 
parameters 

0,n*v(u»v)  “ ««c(|0  v(*)l  : os  2^  p-l}  (1.8) 

max  u,v 

and 

9 (u,v)  *max{|®  (2)J  : 0 5 2 5 p-l)  (1.9) 

max  u,vv 

are  small.  The  same  result  is  obtained  for  q ■ + 1. 

Sot  only  9 (u,v)  itself  but  also  the  number  of  times  |s  (2| 

J max  u,v 

takes  on  this  maximum  value  — when  2 0,1,...,  p-l  --  is  of  Interest. 

Let  jixli  denote  the  cardinality  of  X,  then  we  define 

Lc  - |{{2  : ie  (2)|  » 0 (u,v>)  ; 0 5 & £ p-l}l!  . (1.10) 

u |V  max 


Furthermore  one  defines  for  a set  S of  sequences, 


ft.'.-  *js.y  J5t?v.«-  <‘WUJ 
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0 » max{9  (u,v)  : u€s  , v€s}  . (1.11) 

c max 

A A 

The  parameters  Lc  and  8 are  defined  as  in  (1.10)  and  (1.11)  respectively 

c 

with  0 replaced  by  3. 

Purs  ley  (1976a)  pointed  out  that  the  worst-case  probability  of 

detection  error  P can  be  written  as 

max 

Pmax  “ 1 ’ $(  1 “ (A/pW2£b/N0‘)  (1.12) 

where  $ is  the  standard  Gaussian  cumulative  distribution  function,  S,  * FT 

D 

A 

the  energy  per  data  bit  and  A *»  max{0  ,9  }. 

C 0 

For  the  autocorrelation  functions  we  define  the  parameters 

9 (u)  * max{|9u(Z)!  ,lS|S  p-l]  (1.13) 

La  - II U : |ett(i)|  - e^Ctt)  ;152<  P-i}||  (1.14) 

and  8 m max{9  (u)  : ucS)  (1.15) 

a max  ^ 

A A A 

The  parameters  3^(tr) , La  and  9^  are  defined  as  in  (1.13),  (1.14)  and 

A 

(1.15)  respectively,  with  9 replaced  by  9.  In  addition  we  define 

C^u)  - maxr>ru(X)|  & < p-l)  (1.16) 

and 

C (u,v)  « max{|c  (/)]  : 1-p  ^ ^ p-l)  (1.17) 

max  u,v 

1.4.  Average  performance 

The  average  probability  of  e*.ror  as  well  as  the  average  signal- 


to-noise  radio  at  the  receiver  output  are  important  measures  of  the 
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average  system  performance  given  a set  of  signature  sequences.  The  former 
is  hard  to  compute  but  Yao  (1976)  obtained  upper  and  lower  bounds  on  its 
value.  The  evaluation  of  those  bounds  requires  extensive  knowledge  of 
the  aperiodic  cross-correlation  functions. 

Assume  the  phase  shifts  «p^,  time  delays  and  the  data  symbols 
b^,  q and  b^  1 ^ k ^ K,  k ^ i to  be  mutually  independent  random 
variables.  Also  ass  time  that  cp^  and  are  uniformly  distributed  over 
[0,2Tr]  and  [0,T]  respectively  and  the  data  symbols  b^  ^ take  on  values 
+1  or  -1  with  equal  probability  for  k f i.  By  obtaining  the  mathematical 
expectation  and  the  variance  of  Z^CT) , Purs  ley  (1974)  showed  that  the 
average  (power)  signal-to-noise  ratio  at  the  i-th  user  receiver  output  in 
an  asynchronous  SSMk  system  can  be  expressed  as 


SMa 


(1.18) 


where  the  asynchronous  interference  Q equals 

H 


K 

Q - (6P3)’1  £ t‘(u(1),u(k))  (1.19) 

3 k-1 

k/i 


with  the  interference  parameter  r(u,v)  defined  as 


r(u,v) 


P-1 

E 

/■1-p 


Uc  *(«  + c a) c .d+i)} 


U,v 


u.v 


U.v 


(1.20) 


Yao  showed  that  Pe  * 1 - $ (VSNRa ) is  a very  good  approximation  of  the 
actual  average  error  probability  for  many  practical  values  of  p and  K. 
Clearly,  detailed  knowledge  of  r(u,v)  for  prospective  sets  of  signatures 
is  important  in  the  performance  analysis.  In  some  cases  the  product 


<^^****Wi^^,^^^"w*SroPIWSSy*HlgWiS8raSP^BjppBg«p^j3ggBIBHB>mwa 
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Cu  v^)Cu  V(/+1)  <*oes  not  contr*kute  t0  r(u»v)  as  i-8  shown  in  Appendix  A 
for  Barker  sequences  (Barker?  1953)  of  odd  length  p. 

Important  measures  on  the  autocorrelation  functions  are 


M(u) 


p'-  2 

£ 6,  CO 

4=1  U 


(1.21) 


and 


M(u) 


P-1 

Z 

4*1 


aO 

8Z(A) 

u 


(1.22) 


which  can  be  used  — in  addition  to  the  worst-case  performance  measures  — 
as  sieves  for  sequences  with  good  synchronization  capabilities.  Observe 
that 


M(u)  + M(u)  = 4S(u)  (1.23) 

where 

P-!  2 

S(u)  = E C (4)  , (1.24) 

4=1  u 

Parameter  S(u)  is  called  the  side  lobe  energy  of  a sequence  and  was 
previously  considered  by,  among  others,  Lindner  (1975)  and  Golay  (1977). 


1.5.  Outline  of  the  study 

Chapter  2 investigates  the  various  sequence  parameters  described 
above,  for  random  binary  sequences  for  which  sequence  length  p grows  very 
large.  An  approximation  of  SNR  is  obtained  which  is  very  accurate  for 
typical  values  of  K,  p and  ^,/Nq* 

Chapter  3 discusses  the  code  parameters  for  maximum-. length 
sequences  generated  by  primitive  polynomials.  New  sets  of  sequences  with 
optimal  autocorrelation  as  well  as  cross -correlation  properties  are  obtained. 
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Chapter  4 investigates  the  relationship  between  the  aperiodic 
correlation  functions  of  maximum -length  sequences  and  the  primitive 
polynomials  which  generate  those  sequences.  In  particular,  the  third 
central  moments  of  the  odd  correlation  functions  are  obtained  and  compared 
with  actual  sequence  data. 

Chapter  5 investigates  sets  of  sequences  generated  by  products 
of  primitive  polynomials  such  as  Gold  sequences  (Gold,  1967)  and 
Kasami  sequences  (Kasami,  1966).  A method  --  based  on  Gauss'  products  of 
cyclotomic  cosets  ~ — is  given  which  yields  large  subsets  of  sequences 
with  better  correlation  properties.  The  relevant  correlation  parameters 
of  a large  number  of  good  subsets  are  obtained. 
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CHAPTER  2 

CORRELATION  PARAMETERS  OF  RANDOM  BINARY  SEQUENCES 

This  chapter  considers  the  asymptotic  behavior  of  the  various  code 
parameters  described  in  Chapter  1,  for  random  binary  sequences. 

2.1.  Aperiodic  correlation  parameters 

In  a practical  realization  of  a SSMA  system,  the  signature 
sequence  length  p is  constrained  for  obvious  technical  reasons.  Notwith- 
standing this  fact,  it  is  still  of  interest  to  study  the  asymptotic  behavior 
of  random  binary  sequences  for  which  the  sequence  length  p grows  very  large. 
By  random  binary  sequences  we  mean  binary  sequences  of  independent  identi- 
cally distributed  random  variables  u^,  for  which  Pr{uj=+l}  * Pr(u.«-l}  =*» 

p-l-Z 

Let  0^4^  p-1.  Then,  C (Z)  * £ u v ..  Suppose  |r|  ^ p-Z. 


u,v 


j=0 


j’j+r 


In  order  that  v(Z)  “ r,  it  must  be  that  u^  * v^  f°r  exactly  j(p-Z+r) 
integer  values  of  j in  the  range  0^  j ^ p-l-Z.  Hence,  if  p-Z+r  is  even, 
there  are 

! P-4  \ 

Mi.P.r)  - ^(p./4r)j 

p 2 

choices  for  (Uq,u^, . . . . Since  there  are  lv  choices  for  v and  2 

choices  for  (u  u ,),  there  are  a total  of 

p-*  P“t 

h(Z,p,r)  - 2P  2l  b(X,p,r)  (2.1) 


sequence  pairs  (u,v)  for  which  C (Z)  ■ r,  provided  p-Z+r  is  even.  In 

u,v 

the  special  case  that  u*v,  there  are 
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h'(jJ,p,r)  - 2X  b(A,p,r),  Ui<p-1 


sequence  pairs  (u,u),  i.e.,  sequences  u,  for  which  C u(X)  - r,  l 0. 
Notice  also  that  h(0,p,p)  - 2P. 

Recall  that 

C (u)  ■ max  {|c  (i  ) f : 1 ^ l s p-l} 
max  1 u 1 


and 


C (u,v)  • max  {£  (u,v),  C (v,u)} 

max  max^  * " maxv  ’ J 


where 

“ max  ClC  „<1)|:  0 S l ± p-l}. 

UulX  U|  V 

Moon  and  Moser  (1968)  showed  that  for  e > 0, 


108  Wa> 

h log  p 


* e 


(2.2) 


for  almost  none  of  the  binary  sequences  of  length  p,  i.e.,  for  only  a fraction 
y(p)  of  the  sequences,  where  Xim  y(p)  ■ 0.  Or,  alternatively, 

p «•  ® 


Pr 


log  C (u) 
° max 

k log  p 


\ 


1 


s Y(P) 


(2.3) 


for  u chosen  at  random  from  the  set  of  all  possible  sequences.  In  Appendix 
B we  show,  analogous  to  the  proof  of  (2.2)  by  Moon  and  Moser,  that  if  u and 
v are  drawn  at  random  from  the  set  of  all  2^p  sequence  pairs,  then 


13 


Pr 


where  lim  v(p)  * 0. 

p •*  ® 


log  C (u,v) 
_ maxv  * 

h log  P 


i € 


s v(p) 


) 


(2.4) 


2.2.  Bounds  on  the  set  size  of  'good'  sequence  pairs 

A consequence  of  the  expressions  (2.3)  and  (2.4)  is  the  following. 
Fix  e £ 1.  The  set  of  all  22p  sequence  pairs  is  purged  of  all  pairs  for  which 
the  event  in  (2.4)  occurs.  This  gives  a set  A^(p)  which  has  cardinality 
||A1(p)||.  The  expected  value  of  j|A^ (p)j|  is  lower-bounded  by 

E||A1(p)H  * 22p(l  - Y(p))  . 

Hence,  there  exists  at  least  one  set,  say  A2(p),  with 

||A2(p)H  * 22p(l  - Y(P))  . 

Notice  that  A2(p)  does  not  contain  any  sequence  pair  of  the  form  (u,u) 
because  C (u.u)  ■ p.  This  problem  is  easily  resolved,  however,  because 
expression  (2.3)  implies  that  there  exists  at  least  one  set,  say  B2(p)  with 

l|B2(p)i|  ^ 2P(1  - Y(P». 

The  union  of  A£(p)  and  B2(p)  yields  a set  of  sequence  pairs  with  cardinality 
lower-bounded  by 

i|A2(p)  U B2(p)||  2 22p(l  - ^(p)  + 2’p(l  - Y(p) )) 


(2.5) 


for  which  C (u)  as  well  as  C^^Cu.v)  are  contained  within  the  range 

[p^1"®),  p^1+e:>].  Furthermore,  ||A2(p)  U B2(p)||  grows  exponentially  with  p 
as  22p  because  A^Cp),  which  is  the  compliment  of  A2(p),  and  h_.ice 
||A2(p)  U B2(p)HC#  contains  only  a vanishingly  small  fraction  of  pairs,  i„e., 

2“2pH Ac (p )1J  i y(p)  and  lim  y(p)  ■ 0.  Above  result  shows  that  there  exists  at 

^ P -♦  00 

le*\st  one  very  large  set  of  sequence  pairs  (u,v)  for  which 

£ C (u)  £ , e > 0,  u ■ v 

* maxv 

and 

p^~*  S c (u,v)  £ p**,  e > 0,  u i v . 
r max x 9 

The  size  of  this  set  grows  exponentially  with  p as  22p. 

Discussion 

The  growth  of  the  size  of  at  least  one  set  of  sequences  u,  as  subset 

of  A2(p)  U B2(p),  for  which  then  (2.6)  and  (2.7)  will  hold  simultaneously, 

remains  an  unsolved  problem.  Schneider  and  Orr  (1975)  consider  the  cardinality 
of  a set  A(p)  obtained  by  purging  the  set  of  all  2P  sequences,  of  all  sequences 
which  violate  the  upperbound  of  (2.6)  and  which  form  pairs  which  in  turn 
violate  the  upperbound  of  (2.7).  Their  use  of  the  upperbound  only,  is  here  of 
minor  importance.  The  purging  method  itself,  however,  has  a disastrous  effect 
on  the  lowerbound  of  the  expected  value  of  j|A(p)j|.  Suppose,  for  example,  that 


(2.6) 


(2.7) 
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a certain  sequence  u will  form,  with  many  other  sequences  v,  sequence  pairs 
(u,v)  which  violate  the  (upper)  bound  in  (2.7).  Schneider  and  Orr  — see  their 
equation  (20)  — not  only  remove  sequence  u but  all  the  other  sequences  v as 
well,  from  the  set  of  2P  sequences.  As  a result,  the  lower  bound  of  the 
expected  value  of  ||A(p)||  -fill  not  grow  exponentially  with  p. 


2.3.  Asynchronous  interference  parameters 

Let  u and  v be  two  random  binary  sequences  (not  necessarily  distinct). 
With  |2j  £ p-1  and  |r|  £ p - |ij,  equation  (2.1)  implies 


Pr{c  „(i) 


u.v 


r} 


1 p - Ul 

P ~ Ul+r 
2 


2-<p-UI> . 


(2.8) 


Thi?  probability  mass  function  implies  a moment  generating  function 


«c(t> 


Sfexp[tCu^v(i )] J 

p-2 

n B{  exp[  tu,  V ^33  0 £ l £ p-1 

J-0  J 

p-fc£ 

II  E[exp(tu.  , v 33  1-p  £ l '<  0 

Vj-0  J ‘ 


p-UI 

- n ^(exp(t)  + exp(-c)) 
j-0 


(cosht) 


P*  t 


(2.9) 


Hence  the  first  four  moments  of  C (l)  are 

u,v 
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aM  (C)  n-UI  -1 

E{C  (A)}  - — f- — I - (P“U|)(C08  her'  sin  he.  -0, 
u»v  dC  t-0  't-0 


, 32Mr(t) 

Ef^U))  -TT“ 


| - zf^HWos  ht)p"lil'2sin2ht| 

t-0  v * ' t-0 


at 


+ (p-|x|)(cos  ht)p~^”*cos  ht|  * p-|x|  . 


t-C 


(2.10) 


In  a similar  manner. 


E{C  (X)} 
u,v' 


a Mc(t> 


at- 


t-0 


and 


e{c!  a)} 


s4m  (t> 


u,v 


at"  't-o 


4(P-UI)+(P-W)  + (p.U|, 


3(p-|xl)2  - 2(p-|x| ) . 


(2.11) 


Furthermore,  for  L + n, 


p-l-f£  p-*l-ta 

E{C  (X)C  (o)/X  < 0,m  < O)  - Z Z Efu.v  uv  } • 0. 
u.v'1  u,vv  j-0  n-0  J ***  ^ 


In  fact,  ic  is  easy  to  show  chat 


E{C  (X}C  (o)] 
u,v  u,v  4 


- 0 ¥ X,¥  a,  L + o 


and 


E(c.  (DC  .»}  - 0 


Y I,¥  a,  £ + n 
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while 

" (P"lAl)(P"lmt)  Y l’*  m»  t t m. 

Above  results  enable  us  to  calculate  the  mathematical  expectation  and 
the  variance  of  the  asynchronous  interference  Q as  defined  in  (1.19). 

First  we  obtain  the  first  and  second  ancient  of  the  interference  parameter 
r(u,v). 

P-1  ? 

E{r(u,v)}  2 E{2C^  v(i)  + ^^(^>^^(^+1)} 

P-1 

- 2 2(p-|i|) 

X-l-p 

- 2p2  (2.12) 


and 


E[r2(u,v)}  - E{  ?2  2C2  v(2)}2  +E{  P2  cu  v(X)Cu  va+l>}2 
i»l-p  ’ i-l-p  ’ * 


p-l 


P-1 


+ E{2  2 2C  U)  2 C (o)C  (a+i»  . 
A-l-p  u»  a~l-p  u'  u,v 


(2.13) 


Sow 


P-l 


P-l  •> 


P-l  P-l 


S(  E 2c‘  - 4E(  E c;  «))  +«C  £ E c‘  U)C  (.» 

i-l-p  * Z*l-p  * £-l-p  a-l-p  ' 

It  a 


p-l  o 

-4  2 [3(p-|x|>  - 2(H*I>3 

2“l-p 

r P”1  i5  P-l 

♦ 4 s (p-|i|)i  -4  l (p-U!r 

^tm  i-p  - tm  i-p 


.4  16  3 .1,8 

-4p  -t  “ p - 8p  r j p 


(2.14) 
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and 


P-1 


,2  P;1  .2  ,.„2 


E ( E cu,,«)c„,,«+1>!  - Ef,  ? Cu,v^)Cu,v«+1>J 

i.“l-p  * * i-l-p 


+ Et  ' V V Cu,vW)Cu,v<i+l>Cu,v<”>Cu,v<m+1« 

l»  1-p  m*>l-p  ’ 

i + m 

p-1 

- e (p-UIxp-MI) 

£“l-p 


2 3 2 

3 P ’3P 


(2.15) 


while 


P-1  P-1 

E{2  E 2C,  (i)  E C (m)C  (m+1)}  - 0 . 

A-l-P  U>  m-l-P  ’ ’ 


£*l-p  ' m-i-p 

Hence,  substituting  (2.14),  (2.15),  and  (2.16)  into  (2.13)  gives 


(2.16) 


0 / 0 
E{r  (u,v)}  * 4p  + 6p  - 8p  + 2p  . 


(2.17) 


In  a straightforward  manner  one  obtains  with  (2.12) 


K 


(6P3)-1  Et  £ .<«<».<»)} 


k-1 


(6p3)-1(K-l)  E{r(u,v)3 


(3p)  "(K-1) 


(2.18) 


while 
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var(Qa)  - E(Q 2)  - E2Qa 

■ (36p6)-1(K-l)  var(r(u,v>) 

thus 

var(Qa)  - (36p6)’ V-l.)(6p3  - 8p2  + 2p) 

and  therefore 

var(Qa)  « (6p3)  1(K-i>  (2.i9) 

for  large  sequence  length  p. 

2,4.  Asynchronous  interference  versus  synchronous  Interference 

A system  is  considered  to  be  synchronous  when  the  relative  shift  t 
between  the  signature  sequences  equals  zero.  The  RF-phases  of  the  signal 
carriers,  however,  are  still  assumed  to  be  independent  random  variables, 
uniformly  distributed  over  the  range  [Q,2tt]  . The  signal-to-noise  ratio 
at  the  output  of  a correlation  receiver,  synchronized  (frame  and  bit)  with 
its  own  signature  sequence  equals 

f *Q  "i  -i 

SNRs  • hr + V <2'20> 

b 

where  the  synchronous  Interference  Qg  equals 

Q - (2p2)-1  E r <u(i),u(k))  (2.21) 

3 k-1  8 

k*i 

2 

where  rg(u,v)  ■ v(0).  Wolf  and  Elspas  (1965)  derived  EQg  and  var(Qg) 


for  random  binary  sequences, 
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EQS  - ( 2p) ”1 (K-l ) (2.22) 

and 

var(Qg)  - (2p2)"1(K-l) . (2.23) 

Of  course,  in  a practical  synchronous  system,  one  may  want  to  consider  an 
orthogonal  set  of  signature  sequences  for  which  v(0)  « 0,  as  long  as  the 
synchronization  requirements  can  be  met.  From  (2.18)  and  (2.22)  we  have 
the  well-known  result,  reported  by  Harris  (1973), 

EQa  “ f EQS  . (2.24) 

In  addition,  however,  we  have  now  also 

var(Qa)  - (3p)**1[  1 - J p"1  + J p"2]var(Qg) 

or 

var(Qa)  w (3p)"1var(Qs)  (2.25) 

for  large  sequence  length  p. 

Let  the  asynchronous  fluctuation  ratio  R be  defined  by 

Ra  - 10  log10{[EQa  - V^(Qa)]'1[EQa  + v^7(Qa)]}.  (2.26) 

The  synchronous  fluctuation  ratio  Rg  is  defined  as  in  (2.26)  with 
replaced  by  Q . Substitution  of  EQ  and  var(Q  ) into  (2.26)  gives 

3 a Q 
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Ra  - 10  log10{[K-l-3(K-l)%(6p)"^]"1[K-l+3(K-l)%(6p)^3}  (2.27) 

while  a similar  substitution  with  EQ  and  var(Q  ) gives 

S 8 

Rs  - 10  log10{[K-l-(2(K-l))^]"1[K-l+(2(K-l))^}  (2.28) 

Expressions  (2.27)  and  (2.28)  indicate  a rather  small  value  of  R if 

& 

compared  with  R . Furthermore,  R does  not  depend  on  the  sequence  length  p 
s s 

but  R decreases  steadily  when  p increases.  The  graph  in  Figure  2 gives 

Si 

R and  R as  functions  of  K for  a number  of  sequence  length  p ■ 2n-l.  For 

S Si 

example,  with  K ■ 8 and  p ■ 127,  R * 0.36  dB  while  R “5.1dB*  It  should 

be  noted  that  a fluctuation  of  Qa  does  not  produce  the  same  fluctuation  in 

SNRa.  For  example,  with  10  log^E^/Ng)  *■  10  dB,  a +7var(Qa)  fluctuation 

of  EQa  results  in  an  (approximately)  +0.05dB  fluctuation  of  SNRa  while  a 

+ Vvar(Q  ) fluctuation  of  EQ  results  in  an  (approximately)  + 0.85  dB 

fluctuation  of  SNR  . 

s 

Thus,  in  the  analysis  and  preliminary  design  of  an  asynchronous  SSM& 
system,  the  approximation 

s®a»{^  + ^-1  (2.29) 

is  very  accurate  for  typical  values  of  K,  p and 


CHAPTER  3 


CORRELATION  PARAMETERS  OF  M-SEQUENCES 

In  a practical  spread  spectrum  multiple  access  system,  the  signa 
ture  sequences  will  not  be  selected  at  random  from  the  set  of  all  possible 
sequences  of  a certain  length  p.  The  sequences  must  possess  certain 
qualities,  otherwise  the  SSMA  system  will  not  function  properly.  A 
large  class  of  sequence  with  a number  of  interesting  properties  are  the 
maximum  -length  shift  register  sequences  or  m-sequences.  They  have  been 
studied  extensively  in  the  literature  by  Zierler  (1959)  and  Golomb  (1967) 
and  others. 


3.1.  Introduction  to  m-sequences 

It  is  convenient  to  distinguish  in  our  notation  between  sequence 
elements  u^  6 f -1 , 1}  and  sequence  elements  € {1,0}  which  are  related  by 


Uj  - (-  1)  J • 


(3.1) 


A binary  m-sequence  p,  of  period  p - 2-1,  is  a sequence  which  oatisfies  a 
recurrence  relation  of  the  form 


n 


V> 


fi  h-ta-i’  J"0’1’-*- 


j+u  -i  -J40. 


(3.2) 


where 


f(x)  - fQxn  + flXn_1  + ...  + fn.lX  + fn 


is  a primitive  polynomial  of  degree  n over  GF(2),  the  binary  alphabet 
{0,1}  with  addition  modulo  2.  A polynomial  of  degree  n is  primitive  if  it 


divides  xm  - 1 for  m - 2n  - 1 but  not  for  any  m < 2 -l-  The  roots  of  a 

primitive  polynomial  of  degree  n are  primitive  elements  of  the  extension 
field  GF(2n),  i.e.,  they  have  order  2n  - 1,  and  every  nonzero  element  of 
GF(2n)  can  be  written  as  a power  of  some  primitive  element  P.  We  denote 
the  minimal  polynomial  which  has  Pq  as  a root  by  fq(x).  Peterson  and  Weldon 
(1972)  give  extensive  tables  of  primitive  polynomials  up  to  degree  34. 

Polynomial  f(x)  represents  an  n-stage  linear  feedback  shift 
register  where  £Q  - f n - 1 and  for  0 < i < n,  f£  - 1 if  there  is  a feedback 
tap  connected  to  the  ith  stage  of  the  register  and  £±  * 0 if  not.  An  example 
of  a shift  register  represented  by  f(x)  * x^  + x + 1 is  given  in  Figure  3. 


/^j+4 


FP-5391 


Figure  3.  Shift  register  for  f(x)  - x + x + 1 ; P » 15. 


Let  Tu  denote  the  left  cyclic  shift  of  sequence  u,  i.e. 
TU  - T(Un»^i» * ^2 * »U'p-l,M’0'> 
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Of  course,  Tu  is  equivalent  with  Tu.  Notice  that  equation  (3.2)  implies 

Tn(i  + fx  Tn-1U  + + fn-1  Tpp  + » - 0 (3.3) 

where  0 denotes  the  all-zero  sequence. 

Let  tr(x)  denote  the  trace  of  x,  defined  by 

2 «n-l 

tr(x)  * x + x + + x . (3.4) 

Then  for  each  nonzero  element  ijr  in  GF(2n),  there  exists  a solution  U of 
the  recurrence  relation  (3.2)  specified  by 

pt.  » tr(tpj)  . (3.5) 

A proof  can  be  found  in  Lint  (1973).  The  various  solutions  P- , p.  #0, 
of  (3.2)  are  simply  shifted  versions  of  each  other!  i.e.,  the  sequences 
u,  TP*,  T^u,....,TP  ^p.  are  in  the  same  equivalence  class  with  respect  to 
shift  operator  T;  pi  is  called  the  cycle  representative.  The  m-sequence 
corresponding  to  t a 1 satisfies  the  important  property 

P-j  * ^ (3.6) 

and  is  3aid  to  be  in  its  natural  orientation  or  characteristic  form. 

Henceforth  pi  (or  u)  denotes  the  m-sequence  in  its  natural  orientation, 
k , k 

Furthermore,  T pi  corresponds  with  t a P • For  each  degree  n,  up  to  n * 168, 
Willett  (1976)  has  computed  the  characteristic  form  of  one  m-sequence. 

Let  v be  any  other  m-sequence  of  the  same  length  as  u.  Suppose 
v satisfies  a recurrence  relation  specified  by  f^(x)  and  let  v be  in  its 
natural  orientation,  then 
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’j  ‘ “qj  ’ VJ  • (3-7) 

Equation  (3.7)  implies  that  the  characteristic  m-sequence  v can  readily 
be  determined  from  the  characteristic  m-sequence  u.  Table  1 gives  all 
characteristic  m-sequences  up  to  and  including  length  p a 255.  The 
primitive  polynomials  f(x)  are  denoted  in  the  usual  octal  notation,  for 
example 

f(x)  * x^  + x + 1 = O.x^  + l.x^  + O.x^  + 0.x2  + 1.x  + 1 

3 [0  1 0 0 1 1]  s o 2 3 . 

The  initial  start  position  or  loading  of  the  shift  register  is  also  denoted 
in  octal  notation 


U * (Wq.^.  ) a (0,0,0, 1,0,0,  ) 5 0 4 6 5 . 

Of  course  the  last  two  digits  of  0465  are  here  redundant.  Henceforth,  all 
loading  of  shift  registers .will  be  given  in  octal  notation  except  when 
indicated  otherwise,  by  means  of  an  asterisk,  e.g.,  as  in  Table  4. 


3.2.  The  trinomial  structure  of  m-sequences 

Another  important  property  of  m-sequence  u is  the  so  called 
shift  - and  - add  property,  i.e., 


+ 


j+r  "j+Jt  1 


(3.8) 


for  some  r and  It  0.  In  terms  of  the  corresponding  m-sequence  u over 
{+1  , - 1 } , equation  (3.8)  becomes  uj+r  * uj+iuj  ^or  30me  r anc*  * ^ 


Table  1 
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Characteristic  m-sequences  of  length  p = 15,  31,  63,  127  and  255. 


Poly. 


Loading 


Poly. 


Loading 
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which  will  be  used  extensively  in  Chapter  4.  Equation  (3.8)  can  also 
be  expressed  in  terms  of  the  cyclic  shift  operator  T, 

+ T^tt  + M-  ■ 0 . (3.9) 

Recall  that  m-sequence  M-  already  satisfies  a recurrence  relation  as  in 
(3.2),  specified  by  a minimal  polynomial  f(x).  Hence  (3.9)  specifies 
for  each  l i 0 a binary  trinomial  of  the  form 

xrW+xX+l  (3.10) 

which  should  be  divisible  by  f(x). 

Let  1]  denote  the  cyclotomic  coset  of  integers  x'  mod  p 
containing  x as  smallest  element,  i.e., 

Tlx  * {x‘  : x*  = 2^x  mod  p,  x ^ x'  ; j * 0,1,  . . . . ) . (3.11) 

In  some  cases  T|  will  be  denoted  as  [x]  or  x if  so  indicated.  The 
cyclotomic  cosets  of  integers  modulo  p,  up  to  p = 255  are  tabulated  in 
Table  2.  On  each  horizontal  line  a coset  and  its  reciprocal  are  grouped 
together,  whenever  such  a reciprocal  exists. 

The  trinomials  in  (3.10)  have  algebraic  properties  which  consider- 
ably reduce  the  effort  of  cinding  r(/)  for  each  value  of  i.  Most  important 

is  the  property  that  whenever  l € and  l'  € T]  , than  r(i)  € '0  and 

x x y 

rCi')€T]  , for  some  x and  y.  For  example,  this  can  be  observed  from  Table  3 
where  the  trinomials  xr<^  + %’  + 1,  divisible  by  f(x)  a + x + 1,  are 
tabulated.  A more  detailed  discussion  on  trinomials  can  be  found  in 
Lindholm  (1968).  In  Chapter  4 and  5 we  will  use  the  trinomial  structure  of 
the  m-sequences  extensively. 
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Table  2.  Cyclotomio  oosets  of  integers  mod  p = 15,  31,  63,  127  and  255. 


p = 15 


p = 31 


P = 63 


p = 127 


1 

2 

4 

8 

7 

14 

13 

11 

3 

6 

12 

9 

0 

5 

10 

1 

2 

4 

8 

16 

15 

30 

29 

27 

23 

3 

6 

12 

24 

17 

7 

14 

28 

25 

19 

5 

10 

20 

9 

18 

11 

22 

13 

26 

21 

0 

1 

2 

4 

n 

0 

16 

32 

31 

62 

61 

59 

55 

47 

3 

6 

12 

24 

48 

33 

15 

30 

60 

57 

51 

39 

5 

10 

20 

40 

17 

34 

23 

46 

29 

58 

53 

43 

7 

14 

28 

56 

49 

35 

9 

18 

36 

27 

54 

45 

11 

22 

44 

25 

50 

37 

13 

26 

52 

41 

19 

38 

21 

42 

0 

1 

2 

4 

8 

16 

32 

64 

63 

126 

125 

123 

119 

111 

95 

3 

6 

12 

24 

48 

96 

65 

31 

62 

124 

121 

115 

103 

79 

5 

10 

20 

40 

80 

33 

66 

47 

94 

61 

122 

117 

107 

87 

7 

14 

28 

56 

112 

97 

67 

15 

30 

60 

120 

113 

99 

71 

9 

18 

36 

72 

17 

34 

68 

55 

110 

93 

59 

118 

109 

91 

11 

22 

44 

88 

49 

98 

69 

29 

58 

116 

105 

83 

39 

78 

13 

26 

52 

104 

81 

35 

70 

23 

46 

92 

57 

114 

101 

75 

19 

38 

76 

25 

50 

100 

73 

27 

54 

108 

89 

51 

102 

77 

21 

42 

84 

41 

82 

37 

74 

43 

86 

45 

90 

53 

106 

85 

p = 255  0 


1 

2 

4 

8 

16 

32 

64 

128 

127 

254 

253 

251 

247 

239 

223 

191 

3 

6 

12 

24 

48 

96 

192 

129 

63 

126 

252 

:49 

243 

231 

207 

159 

5 

10 

20 

40 

80 

160 

65 

130 

95 

190 

125 

250 

245 

235 

215 

175 

7 

14 

28 

56 

112 

224 

193 

131 

31 

62 

124 

248 

241 

227 

199 

143 

9 

18 

36 

72 

144 

33 

66 

132 

111 

222 

189 

123 

246 

237 

219 

183 

11 

22 

44 

88 

176 

97 

194 

133 

61 

122 

244 

233 

211 

167 

79 

158 

13 

26 

52 

104 

208 

161 

67 

134 

47 

94 

188 

121 

242 

229 

203 

151 

15 

30 

60 

120 

240 

225 

195 

135 

17 

34 

68 

136 

119 

238 

221 

187 

206 

19 

38 

76 

152 

49 

98 

196 

137 

59 

118 

236 

217 

179 

103 

157 

21 

42 

84 

168 

81 

162 

69 

138 

87 

174 

93 

186 

117 

234 

213 

171 

23 

46 

92 

184 

113 

226 

197 

139 

29 

58 

116 

232 

209 

163 

71 

142 

25 

50 

100 

200 

145 

35 

70 

140 

55 

110 

220 

185 

115 

230 

205 

155 

27 

54 

108 

216 

177 

99 

198 

141 

39 

78 

156 

57 

114 

228 

201 

147 

37 

74 

148 

41 

82 

164 

73 

146 

91 

182 

109 

218 

181 

107 

214 

173 

43 

86 

172 

8 9 

178 

101 

202 

149 

53 

106 

212 

169 

83 

166 

77 

154 

45 

90 

180 

105 

210 

165 

75 

150 

51 

102 

204 

153 
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Table  3.  Trinomials  x ^ ' + x +1  divisible  by  f(x)  ■ xH  + x + 1 


X^  + X + 1 

13  , 6 , , 

12  . 11  , , 
X + X + 1 

8 2 . 
x + x +1 

X + x +1 

11  . 12 

X + X + 1 

14  . 3 . _ 

X + X + 1 

X2  + X8  + 1 

6 . 13  . , 

X + X +1 

X + x^  + 1 

X7  + X9  + 1 

3 . 14  . , 

X +x  +1 

X + X + 1 

X5  + X10  + 1 

3.3.  Autocorrelation  functions 

The  autocorrelation  function  of  a signature  sequence  plays  a key 
role  in  obtaining  word  synchronization  in  the  correlation  receiver  and 
in  reducing  the  effects  of  multipath  interference.  The  periodic  auto- 
correlation function  8U(T)  " “ 1 for  all  f ^ 0 mod  p for  m-sequences  is 
nearly  ideal  in  this  respect  and  one  of  the  main  advantages  of  m-sequences. 
Of  course  ease  of  generation  is  another  advantage. 

While  thus  8^  (T  u)  * 1,  Vk  for  m-sequences,  the  companion  param- 

* 

eter  9 (T'  u)  defined  in  (1.13)  is  very  sensitive  to  the  selected  cyclic 

max 

shift  of  u.  To  find  cyclic  shifts  for  which  9 (Tku)  is  minimal, 

max 

a computer  search  is  required.  In  the  process  of  searching  the  best  T*, 

a number  of  ties  may  occur  (i.e.,  a number  of  different  values  of  k in 
k n k 

T result  in  the  same  minima  1 9 (T  u)).  Hence,  a second  condition  is 

max 

applied  to  the  already  selected  values  of  k reducing  the  number  of  ties 
considerably.  The  second  condition  is  stated  in  the  following  definition 
which  is  due  to  Massey  and  Uhran  (1969). 
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Definition  1.  An  m-sequence  u'  ■ x u is  auto-optimal  (AO)  — with  respect 

a 

to  a (T  u)  — when  the  following  conditions  are  satisfied  in  successive 
niax 

order: 


a) 

b) 


e(u')  ^ § (A),  Vk 

max  max 


The  cardinality  La  of  the  set 


C^:|l  (A)|  » 0 (u*);  0 < l<  p-1] 

u max 

is  smallest  for  sequence  u* . 


Observation:  It  follows  immediately  from  this  definition  that  an  auto- 
optimal  m-sequence  generated  by  primitive  polynomial  f(x)  of  degree  n has 
a reciprocal  — generated  by  f'(x)  e xnf(l/x)  --  which  is  AO  too. 

Table  4 specifies  for  each  primitive  polynomial  of  degree  n « 7, 
the  loading  or  start  position  Uq . . . . of  the  shift  register  such 
that  the  generated  m-sequence  is  AO.  There  are  two  distinct  cyclic  shifts 

7 

of  the  m-sequence  generated  by  £(x)  * x + x + 1 (203),  both  which  are  AO. 

They  are  indicated  as  203a  and  203b.  Then,  the  reciprocal  polynomial 
7 6 

£(x)  “x  + x +1  (301)  will  generate  two  distinct  AO  ta-sequences  too. 

a 

The  resulting  values  of  9 (u)  and  La  for  u*  are  as  indicated. 

max 

Clearly,  there  is  not  one  unique  set  of  eighteen  ^uto-optioal 
ta-sequences  of  length  p * 127  as  was  reportad  by  Massey  and  Uhran  (1969). 
Table  4 shows  that  a total  of  four  dlstinet  sets  of  eighteen  o-sequences 
each  can  be  selected  by  choosing  two  a-sequences,  which  arc  not  cyclic 
shifts  of  each  other,  from  the  set  (203a,  203b,  301a,  301b}.  Furthermore, 
the  auto-optlsal  a-sequences  of  Massey  and  Uhran  do  not  seeo  to  be 
reciprocal  pairs. 
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Table  4.  Auto-optimal  m-sequencea  of  length  p = 127. 


Poly. 

Loading* 

Poly. 

Loading* 

0mfl„(u) 

max 

A 

La 

S(u) 

211 

0010000 

221 

1001101 

17 

6 

2183 

217 

0000101 

361 

1111111 

15 

12 

2015 

235 

OOCIIOO 

271 

1000101 

17 

10 

2283 

247 

0010111 

345 

0110001 

17 

8 

2255 

277 

1110001 

375 

0101010 

IS 

4 

2295 

357 

1110010 

367 

0110101 

17 

4 

2563 

323 

1110111 

313 

1000111 

17 

4 

2203 

203a 

1101101 

301a 

0010010 

17 

4 

2087 

203b 

0000001 

301b 

1111111 

17 

4 

2403 

325 

0000101 

253 

1101100 

19 

6 

2483 

The  sidelobe  energy  parameter  S(u),  as  defined  in  (1.24)  can  also 
be  used  as  a sieve  for  m-sequences . In  particular,  we  might  use  S(u)  to 
further  distinguish  between  AO  m-sequences,  because  for  m-sequences, 


4 S(u)  - p + 1 


(3.12) 


one  of  the  important  parameters  in  the  direct  sequence  SSMA  system. 

Hence,  it  is  convenient  to  extend  the  definition  of  auto -optimality 
further  by  including  3(u)  as  a next  sieve. 

k* 

Definition  2.  An  m-sequence  U * T u is  indicated  as  AO/LSE  whenever  U 
is  auto-cptimal  and  has  lowest  sidelobe  energy  S(U)  among  all  auto-optimal 
shifts  of  u. 

When  the  m-sequences  generated  by  polynomials  203b  and  301b  are 
deleted  from  Table  4,  the  table  will  give  the  AO/LSE  m-sequences  of  length 


p =*  127. 


33 


3.4.  Periodic  cross -correlation  of  m- sequences ; Golomb's  theorem 

Let  u and  v be  characteristic  m-sequences  generated  by  the 
primitive  polynomials  f^(x)  and  f^(x)  respectively.  Both  m-sequences 
are  constant  over  cyclotomic  cosets,  i.e., 


Uj  - u2j  = xO^)  for  j € T\± 


(3.13) 


and 

where 


TJ  " T2j  ' X<V  J€\t 


X(H 


qi 


(3.14) 


(3.15) 


The  periodic  cross -correlation  0 (t)  is  also  constant  over  eye lo comic 

U,  V 

cosets,  i.e.. 


0 (t‘)  - e (2r')  - 6 ,<1L),  t’€T1  . 

U,vx  U,V  U,V  T T 


(3.16) 


Gold  and  Kopitzke  (1963)  have  computed  the  periodic  cross-correlation 
0 V(T)  f°r  m-sequences  of  length  p ^ 8191.  The  same  data  can  be  obtained 
by  means  of  an  elegant  theorem  of  Golomb  (1968),  which  is  derived  as  follows. 
The  Gauss'  product  of  cyclotomic  cosets  and  ^ is  defined  as 


VV * 0104  p : s ^ \ » t € • 


(3.17) 


Hence,  one  can  write  the  following  equality 


E £ (-1) 

ten 


trtf**) 


ren  n 

x y 


(-1) 


tr(8r) 


(3.18) 


Let  jlH  II  denote  the  cardinality  of  cyclotomic  coset  T)  . Then,  from  the 


definition  of  9 (t) , equations  (3.15),  (3.16)  and  (3.18),  one  obtains 

u,v 


Golomb's 
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Theorem  0u,v(V  " jj^jj  2 Xu(T1-A)X(\i)  (3‘19) 

T 

where  X(\\)  « E (-l)**^  (3.20) 

u x y r€T|  T1 
x y 

and  the  sum  in  (3.19)  is  over  all  the  integers  i which  represent  distinct 
cyclotomic  cosets  T|  . Expression  (3.19)  represents  a multiplication  of  a 
vector 

Xq  - (X(Hql))i  (3.21) 

of  'coset  assignments'  to  m-sequence  v,  and  the  normalized  matrix 


X = [ 
~u 


— X (T1  H ))  . 

IIHj  U Xy 


(3.22) 


In  order  to  obtain  X for  some  m-sequence  u of  period  p one 

<“'-U 

needs  the  array  of  Gauss'  cyclotomic  coset  products  £T[  T)  } . Golomb  (1963) 

x y 

specified  (T|^T|  } for  p * 15  — here  reproduced  in  Table  5 --  and  p = 31; 

we  present  {T1  T|  } for  p ® 63  in  Table  6,  where  T|  is  indicated  by  [x] . 
x y x 


An  example  of  Golomb* s theorem 

Let  u and  v be  m-sequences  of  length  p “ 15  generated  by  the 
4 4 3 

polynomials  f.(x)  ■ x + x + 1 (023)  and  f-,(x)  ■ x + x +1  respectively. 

4 * 

The  cyclotomic  cosets  are  “Hq  c (0) , ® {l,2,4,8},  ^ “ (3,6,12,9), 

Tl5  • (5,10}  and  ^ • (7,14,13,11).  The  vectors  of  coset  assignments  to 
m-sequences  u and  v are 


Xi  » (xoy,  x(y.  X(H3),  xoy,  X(T|?)) 
- (l,  l,  -l,  i,  -l) 


\ 

\ 

11,  +2lL+2Tlr 
13  3 

11  +211. +211, 

1 5 7 

4^+^+^+^  . 

s 

\ 

H1+2115+2T7 

4*n0  +3T|3 

\*h 

2H1  +2T15  +H7 

s 

s 

V^7 

V*7 

2”o  +"5 

V^3 

4H0+Hi+H3+tl7 

'2Tli  +21)5  +117 

V^3 

2*n3  +2115  +ti7 

and 

X7  - (xC'Hq),  x(T17),  x(T|3),  xcn5),  x(\)> 
- (i,  -i.  -i,  i»  i)  • 


Evaluation  of  the  Gauss'  products  {11  T)  } in  Table  5 with  (3.20)  and  (3.22) 

x y 

results  in  the  matrix 


5u 


4 

0 

0 

-4 

0 


-4  2 

0 -2 
-2  0 
0 2 
2 0 


Multiplication  of  and  X^  results  in  a vector 


(3 


l,  „0U  ■ (*1.  -1.  3,  7,  -5) 
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38 


or  alternatively 


0 

u,v 


(T) 

V 


-1  if  t€T|0 
-5  if  t€  T) 

3 if  t€T13 
7 if  t€  T)5 
-1  if  t6T17  . 


(3.24) 


Whenever  two  m-sequences,  say  w and  z are  decimations  of 
the  m-sequences  u and  v respectively  such  that  = ug  ^ and  z^  * v , then 


0 

w,z 


(T) 


P-1 

USj  VSj4-ST  " 6U,V(ST) 


(3.25) 


Hence,  to  obtain  all  periodic  cross -correlation  values  for  the  m-sequences 
of  a certain  length  p,  it  suffices  to  determine  one  matrix  X^.  For 
p * 15  and  p * 31,  a matrix  X is  given  in  (3.23)  and  in  Golornb  (1968) 
respectively.  Table  7 and  8 specify  a matrix  X^  obtained  by  evaluating 
the  Gauss'  products  of  cyclotomic  cosets  for  the  indicated  m-sequences 
of  length  p * 63  and  p = 127. 

We  will  encounter  Gauss'  product  of  cyclotomic  cosets  in  com- 
bination with  Golomb's  theorem  again  in  Chapter  5 where  the  periodic  cross  - 
correlation  of  Gold  sequences  and  Kasami  sequences  are  discussed. 


3.5.  Maximal  connected  sets 

Gold  and  Kopitzke  (1965)  obtained  from  their  data,  sets  of 
m-sequences  of  period  p,  called  maximal  connected  sets,  which  are  the 
largest  possible  subsets  of  m-sequences  for  which  any  two  sequences  in 
the  same  set  have  a preferred  three-valued  periodic  cross-correlation 


Table  7.  Gauss*  products  of  cyclotomic  cosets  evaluated  for 
m-sequence  u with  octal  polynomial  103;  p = 63. 


0 1 

1 6 

1 2 

1 2 

-1  2 

1 -2 

1 -2 

-1  2 

1 -2 

-1  -2 

-1  -6 

-1  -2 

1 -2 

-1  2 


3 5 

6 -6 
2 2 
2 -2 

-2  -4 

-2  2 
-2  2 
-2  0 

-2  -2 

2 0 
6 0 
2 0 

-2  -2 

-2  4 


7 9 

6 3 

-2  -1 

-2  -1 
2 1 
2 -1 
-2  3 

-2  -3 
2 -1 
-2  1 
6 -3 
-2  1 
-2  3 

2 1 


11  13 

-6  6 
2 -2 
-2  -2 
0 -2 
-2  2 
-6  -2 
0 2 
2 2 
0 2 
0 -6 
0 2 
6 -2 
0 -2 


15  21 

-6  -2 
-2  -2 

2 2 
0 0 
-2  2 
2 -2 
0 0 
2 -2 
4 0 

0 0 
-4  0 

-2  2 
0 0 


23  27 

-6  3 

-2  -1 
2 -1 
0 -1 
-2  -1 
2 3 

0 3 

2 -1 
-4  -1 
0 3 

4 -1 
-2  3 

0 -1 


31 


-6 


2 

-2 

4 

2 

2 

0 

-2 

0 

0 

0 

-2 

-4 


Table  8.  Gauss'  products  of  cyclotomic  cosets  evaluated  for 
m-sequence  u with  octal  polynomial  211;  p a 127. 


0 

1 

3 

5 

7 

9 

11 

13 

15 

19 

21 

23 

27 

29 

31 

43 

47 

55 

63 

-1 

7 

7 

7 

-7 

7 

-7 

7 

-7 

-7 

-7 

-7 

-7 

7 

-7 

7 

7 

-7 

7 

1 

7 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

7 

-1 

-1 

-1 

-1 

1 

-1 

-1 

3 

-1 

-1 

-1 

-1 

3 

-5 

3 

3 

-1 

3 

-1 

-5 

-1 

-1 

3 

1 

-1 

3 

-1 

3 

-5 

-5 

-1 

3 

3 

-1 

-1 

-1 

-1 

-1 

3 

-1 

-1 

3 

-1 

-1 

-1 

3 

1 

-1 

1 

-1 

-3 

5 

-3 

-3 

1 

3 

1 

-5 

-1 

1 

3 

1 

-1 

1 

-5 

-1 

3 

3 

-1 

3 

-1 

-5 

-1 

3 

-1 

-1 

-1 

3 

-1 

3 

-1 

-1 

-1 

-5 

1 

3 

-3 

-1 

-3 

1 

5 

1 

-3 

-1 

1 

-1 

3 

1 

3 

1 

-1 

-1 

-1 

-1 

-1 

-1 

7 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

-1 

7 

-1 

-1 

-1 

3 

3 

-3 

3 

-3 

-1 

1 

1 

-3 

5 

1 

-5 

1 

-1 

-1 

1 

-1 

-1 

-1 

-5 

3 

5 

-1 

1 

-1 

1 

-3 

-3 

1 

-3 

-1 

1 

3 

3 

1 

-1 

-1 

-1 

3 

-1 

-3 

-5 

5 

-1 

-3 

-3 

1 

1 

1 

-1 

1 

3 

-1 

1 

3 

-1 

-1 

3 

-1 

-3 

-1 

! 

-1 

5 

1 

1 

-3 

-3 

3 

1 

-1 

3 

1 

-5 

-1 

-1 

-1 

-1 

1 

3 

-3 

-1 

1 

-3 

1 

-3 

5 

3 

1 

3 

-5 

1 

-1 

1 

-1 

3 

-1 

3 

-1 

-1 

-1 

-5 

-1 

-1 

3 

3 

3 

-1 

-1 

3 

-1 

-5 

-1 

7 

-1 

-1 

1 

-1 

1 

-1 

1 

1 

1 

1 

1 

-1 

-7 

-1 

-1 

1 

-1 

1 

-1 

-5 

3 

-5 

-1 

-1 

-1 

-1 

3 

3 

-1 

3 

-1 

-1 

3 

3 

-1 

-1 

1 

-1 

-1 

-1 

-1 

3 

3 

-1 

-1 

3 

-1 

3 

-5 

3 

-1 

3 

-5 

-1 

-1 

-1 

-1 

-1 

-1 

1 

-1 

1 

7 

1 

1 

1 

1 

1 

-1 

1 

-1 

-1 

-7 

-1 

1 

-1 

3 

3 

3 

3 

3 

-1 

-1 

-1 

3 

-5 

-1 

-5 

-1 

-1 

-1 

-1 

-1 
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function,  i.e.,  0 (t)  takes  on  the  values  2 L(n+2)/2j 

u ,v 

and  -1  only.  Any  pair  of  primitive  polynomials  which  yields  a sequence 
pair  with  a preferred  three-valued  periodic  cross-correlation  is  called 
a preferred  pair  of  polynomials  as  discussed  in  Pursley  (1976).  Clearly, 
for  preferred  pairs  of  polynomials,  which  generate  m-sequences  u and  v, 
the  peak  parameter  defined  in  (1.8)  equals 

WU>V)  ’ 2L<n+2)/2J  + l • 

For  example,  the  first  six  polynomials  in  Table  4 form  a maximal 
connected  set  S of  six  m-sequenees  of  the  same  length  p * 127  for  which 

®c  * maxlQ^O^v)  : u€s,  v€s}  **  1 + 2 L(6-K2 ) /2j  _ ^ . 

Other  m-sequences  of  length  p * 127,  not  maximal  connected,  can  form 

pairs  with  correlation  peak  values  as  high  as  ®max(u>v)  = 41*  The  size 

of  a maximal  connected  set  of  m-sequences  is  rather  small.  In  fact 

all  sets  have  a cardinality  less  or  equal  to  six  for  p ^ 4095  while  in 

some  cases,  such  as  for  n « 0 mod  4,  the  cardinality  is  zero. 

In  a first  analytical  result,  Gold  (1967)  established  a bound 

on  9 (u,v)  as  a function  of  u and  v. 

max 

Theorem  2:  Let  f^(x)  be  a primitive  polynomial  of  degree  n and  let  P 
be  a root  of  f^(x)  in  GF(2n).  if  ft(x)  is  the  minimal  polynomial  of 
where  t -*  2^-(n+2)/2j  + ^ then  ®max(>u»v)  ^ ^or  sa<lu®nc®s  generated 
by  f1(x)  and  ffc(x). 

Notice  that  ft(x)  does  not  have  to  be  primitive.  Theorem  2 is 
a special  case  of  more  general  results  obtained  by  Kasami  (1966)  for  a 
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large  number  of  values  t * 2 +1.  One  such  a result  is 


Theorem  3:  Let  f^(x)  as  in  Theorem  2 and  ffc(x)  is  the  minimal  polynomial 

of  ft(x)  where  t = 2 +1  and  0 < h < n.  Let  c denote  the  greatest 

common  divisor  of  the  integers  n and  h.  If  n/c  is  odd  then  0y  v(r)  * - 1 

for  2n  -2n  C-  1 values  of  T,  0^  v(t)  55  - 1 - 2^n+c^  for 

2a-c-l.  2 (n-c-2)/2  w1bm  o£  T a’d  9^(t)  . . j + for 

+ 2^n  c values  of  T. 

A more  complete  discussion  is  given  by  Pursley  (1976).  The 
results  of  Gold  and  Kasami  form  a basis  for  the  construction  of  large  sets 
of  sequences  with  good  periodic  correlation  properties  and  are  discussed 
in  Chapter  5. 


3.6.  Aperiodic  correlation  functions 

From  the  analysis  of  the  direct  sequence  SSMA  system  as  pre- 
sented in  Chapter  1,  it  is  clear  that  the  odd  cross -correlation  function 
v(A),  as  defined  in  (1.6)  is  as  important  as  0^  y(4).  As  shown  in 
sections  3.4  and  3.5,  0^  y(4)  displays  certain  regularities  as  a function 
of  j i and  the  polynomials  f(x)  which  generate  u and  v.  The  odd  cross- 

A 

correlation  0^  v(4),  however,  seems  to  be  refrained  from  any  such 
regularities  — and  must  be  computed  for  each  l — hence  the  peak  param- 
eter  9max(u,v) , as  parameter  (u)»  *s  ver7  sensitive  to  the  cyclic  shifts 
of  the  sequences  u and  v. 

Whenever  the  phase  shifts  of  the  sequences  are  already  fixed  due 
to  other  requirements  such  as  a low  value  of  the  autocorrelation  parameter 

A 

®max^’  t*16  cr08S -correlation  parameters  can  readily  be  computed. 
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For  example,  in  Table  9 the  parameters  0 (u),  6 (u,v)  and 

9fflax(u»v)  are  tabulated  on,  below  and  above  the  diagonal  respectively,  for 

the  auto-optimal  m-sequences  of  length  p ■ 127.  These  sequences  are 

indicated  by  the  polynomials  f(x)  which  generate  them.  The  symmetries  in 

the  table  are  due  to  the  fact  that  § (u*,v')  * 9 (u,v)  and  0 (u',v') 

max  max  max 

3 9max^U,V^  w*ienever  u'  an<*  v'  are  reciprocals  of  u and  v respectively 

tc  A 

(see  Section  3.3).  Whenever  v = T u,  6 (u,v)  * 0 (u,v)  « p.  Hence, 

v 7 max  max  r ’ 

those  values  are  omitted  from  the  table.  Notice  also  that  for  the  maximal 
connected  set  S formed  by  the  first  six  m-sequences  in  Table  9, 

0 * max{9  (u,v)  : u € S , v 6 S)  =■  33  . 

c max 

The  interference  parameter  r(u,v)  as  defined  in  (1.20)  for  the 
auto-optimal  m-sequences  of  length  p » 127  is  tabulated  in  Table  10. 

The  auto-optimal  m-sequences  in  Table  10  are,  as  in  Table  9,  indicated 
by  the  polynomials  f(x)  which  generate  them.  Only  one-fourth  of  the 
interference  parameter  values  are  tabulated  because  r(u,v)  ■ r(u,v')  when- 
ever v and  v'  are  reciprocal  sequences.  Notice  the  higher  interference 
between  an  m-sequence  and  itself  or  its  reciprocal.  This  is  a good 
argument,  not  to  use  an  m-sequence  and  its  reciprocal  in  the  same  set  of 
signature  sequences. 

It  is  not  surprising  that  m-sequences,  which  have  many  charac- 
teristics in  cotanon  with  random  binary  sequences,  yield  values  of  r(u,v) 

2 

which  are  close  to  Er(u,v)  - 2p  as  in  (2.12).  In  fact,  all  values 
of  r(u,v)  in  Table  10  --  the  values  on  the  diagonal  excluded  --  are  within 


the  range  [Er(u,v)  + J var(r(u,v))]  for  random  binary  sequences. 
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Table 


9.  Correlation  values  in  the  auto-optimal  set  of  m-sequences;  p = 


2 

2 

2 

2 

2 

3 

3 

2 

2 

3 

2 

3 

2 

3 

3 

3 

3 

3 

3 

2 

1 

1 

3 

4 

7 

5 

2 

0 

0 

2 

2 

6 

7 

4 

7 

6 

1 

0 

0 

5 

1 

7 

5 

7 

7 

» 

3 

3 

3 

5 

1 

1 

1 

5 

5 

7 

3 

1 

1 

3 

u 

b 

a 

b 

211 

17 

33 

31 

27 

27 

31 

39 

27 

35 

25 

37 

45 

33 

47 

33 

35 

29 

37 

29 

27 

217 

17 

15 

31 

31 

27 

29 

31 

37 

35 

31 

45 

35 

23 

37 

33 

29 

29 

31 

43 

29 

235 

17 

17 

17 

27 

31 

29 

29 

35 

27 

29 

33 

23 

37 

29 

33 

27 

31 

33 

33 

31 

247 

17 

17 

17 

17 

31 

29 

33 

29 

31 

33 

47 

37 

29 

35 

33 

29 

29 

23 

25 

53 

277 

17 

17 

17 

17 

19 

27 

33 

33 

25 

39 

33 

33 

33 

33 

37 

31 

29 

29 

29 

35 

357 

17 

17 

17 

17 

17 

17 

31 

39 

29 

27 

35 

29 

27 

29 

31 

41 

31 

31 

25 

35 

323 

17 

17 

41 

17 

17 

17 

17 

41 

31 

41 

29 

29 

31 

29 

29 

31 

37 

27 

27 

27 

203a 

17 

17 

41 

17 

17 

41 

17 

17 

37 

37 

31 

33 

23 

29 

31 

27 

37 

25 

25 

203b 

17 

17 

41 

17 

17 

41 

17 

17 

35 

29 

43 

33 

25 

29 

25 

27 

25 

43 

33 

325 

17 

17 

17 

17 

41 

17 

41 

41 

41 

19 

27 

29 

31 

53 

35 

35 

27 

25 

33 

37 

221 

21 

41 

41 

41 

41 

41 

41 

17 

17 

17 

17 

33 

31 

27 

27 

31 

39 

27 

35 

25 

361 

41 

21 

17 

41 

41 

41 

41 

41 

41 

17 

17 

15 

31 

31 

27 

29 

31 

37 

35 

31 

271 

41 

17 

21 

41 

17 

41 

17 

17 

17 

41 

17 

17 

17 

27 

31 

29 

29 

35 

27 

29 

345 

41 

41 

41 

21 

41 

41 

17 

17 

17 

41 

17 

17 

17 

17 

31 

29 

33 

29 

31 

33 

375 

41 

41 

17 

41 

21 

17 

41 

41 

41 

17 

17 

17 

17 

17 

19 

27 

33 

33 

25 

39 

367 

41 

41 

41 

41 

17 

21 

17 

17 

17 

41 

17 

17 

17 

17 

17 

17 

31 

39 

29 

27 

313 

41 

41 

17 

17 

41 

17 

21 

41 

41 

17 

17 

17 

41 

17 

17 

17 

17 

41 

31 

41 

301a 

17 

41 

17 

17 

41 

17 

41 

21 

21 

17 

17 

17 

41 

17 

17 

41 

17 

17 

37 

301b 

17 

41 

17 

17 

41 

17 

41 

21 

21 

17 

17 

17 

41 

17 

17 

41 

17 

17 

35 

253 

17 

17 

41 

41 

17 

41 

17 

17 

17 

21 

17 

17 

17 

17 

41 

17 

41 

41 

41 

19 

Table  10.  Interference  parameter  r(u,v)  for  the  auto-optimal 
set  of  m-sequences  of  length  p = 127. 


211  217  235  247  277  357 


323  203a  203b  325 


211 

41214  33622  32722  32022  29070 

217 

40222  30954  31406  31446 

235 

42046  32570  33714 

247 

40326  34054 

277 

41294 

357 

323 

203a 

203b 

325 

34394  32978  31486  33114  30250 
29554  32666  34118  35186  32066 
32974  322H6  33106  30766  33310 
32002  29546  30486  34194  33634 
31146  31938  33934  31290  34202 
41958  31258  30806  30626  30770 
39870  33538  32670  33598 
39902  29698  32442 
42358  33830 
42894 


127. 
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In  Appendix  C,  the  correlation  parameters  § (u),  9 (u,v) 

max  max 

and  6 (u,v),  as  well  as  interference  parameter  r(u,v)  are  given  for 

max 

the  AO/LSE  m-sequences  (see  Definition  2)  of  length  p = 31,  63,  127  and 
255. 

Whenever  it  is  possible  to  relax  the  autocorrelation  require- 

A 

ments,  one  might  try  to  minimize  the  cross -correlation  parameter  0 (u,v). 

max 

Even  for  small  sets  of  ~aquences,  however,  the  amount  of  computation 
required  to  find  the  cyclic  shifts  which,  for  example,  minimize  0 

c 

becomes  rapidly  prohibitive.  A discussion  about  the  computational  com- 
plexity of  the  correlation  problem  is  given  by  Pursley  and  Sarwate  (1976). 
Sywyk  (1975)  obtained  some  results  for  auto-optimal  m-sequences  of  length 
p “ 63.  A sub-optimal  result  for  the  above  mentioned  imr*imal  connected 
set  S of  six  m-sequences  of  length  p =°  127  is  reported  in  Table  11.  With 

A A 

the  indicated  binary  loading,  one  obtains  6 * 29  compared  with  9 * 33 

c c 

A 

for  the  auto-optimal  loading.  Notice,  however,  that  9 has  increased 

A A 

from  a previous  low  value  of  8 * 19  (sea  Table  4)  to  8 *23.  The 

3 3 

resulting  values  of  r(u,v)  are  also  given  in  Table  11. 

Finally,  to  obtain  an  indication  of  how  small  8 could  be,  a 

c 

triple  of  polynomials  (211,217,235)  was  selected  from  the  maximal  connected 

A 

set  S.  We  established  that  a lowest  value  8 ° 23  is  achieved  for  two 

c 

sets  of  register  loadings,  here  reported  in  Table  12.  Notice  that  6^  “ 29 

and  8 “21  respectively, 

c 


Table  1 1 . Maximal  connected  set  of  m-3equences  for  which  §c  ■ 29 
and  interference  parameters  r(u,v)  for  this  set. 


Poly. 

Loading* 

§ a <tt) 
max 

217  235 

247 

277  357 

211 

1101010 

19 

217 

1101111 

2? 

211 

32770  28998 

33202 

32870  32290 

235 

1101000 

19 

217 

32810  33014 

33050  32038 

247 

1111111 

23 

235 

30298 

31550  32730 

277 

0111111 

23 

247 

32458  33350 

357 

0101110 

23 

277 

33482 

Table  12.  Optimal  triples  of  m-sequences  { 211,  217,  235}  with  9 

c 

and  interference  parameters  r(u,v)  for  these  triples. 


Poly. 

Loading* 

& (u) 

max 

217  235 

211 

0100000 

21 

217 

1010011 

29 

211 

30494  28733 

235 

1111101 

23 

217 

33362 

211 

0100000 

21 

217 

0010001 

19 

211 

29714  28790 

235 

1110111 

19 

217 

32442 
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CHAPTER  4 

ON  THE  MOMENTS  OF  THE  APERIODIC  CORRELATION  FUNCTIONS 

While  the  relationship  between  the  periodic  correlation  parameters  of 
m-sequences  and  the  polynomials  f(x)  is  rather  well-understood  (see  Chapter 
3),  it  is  less  clear  how  those  polynomials  relate  to  the  aperiodic  correla- 
tion parameters.  In  this  chapter  we  report  our  work  on  one  particular 
relationship  which  had  the  potential  of  being  a possible  sieve  for  m- 
sequences  with  good  aperiodic  correlation  parameters. 

4.1.  Third  moment  problem 

Our  investigations  were  inspired  by  the  work  of  Lindholm  (1968)  who 
was  interested  in  the  weight  distribution  of  M-tuples  of  long  m-sequences. 
He  actually  established  a relationship  between  the  polynomial  f(x)  which 
generates  the  m-sequence  and  the  moments  of  the  M-tuple  weight  distribu- 
tion. Related  results  were  obtained  by  Wainbeig  and  Wolf  (1970)  who 
derived  the  first  six  moments  of  the  M-tuple  weight  distribution  while 
Weathers,  et  al,,  (1974)  obtained  expressions  for  hybrid  m-sequences. 

In  this  chapter  we  consider  the  first  few  moments  of  the  aperiodic 
correlation  function  values  when  the  correlated  m-sequences  have  an  at 
random  selected  cyclic  phase  shift.  In  particular,  we  obtain  the  first 
few  moments  of  the  odd  correlation  function  values.  For  a number  of 
m=- sequences , the  moments  were  evaluated,  and  compared  with  actual  data 
acquired  from  those  m-sequences. 


“V  * 


4.2.  Moments  of  the  aperiodic  correlation  functions 


Let  u * Tu  and  v * Tyv  be  two  cyclic  shifts  of  m-sequences  u and  v 
respectively.  When  we  assume  x and  y to  be  uniformly  distributed  over  the 
integer  values  in  the  range  [0, p-1],  the  moments  of  the  aperiodic  auto- 
correlation function  C~(2)  and  the  aperiodic  cross -correlation  function 

C~  ~(2)  can  be  calculated  in  a relatively  straightforward  manner,  as  is 
U,v 

shown  in  Appendix  D. 

In  Appendix  D it  is  also  shown  that  the  mathematical  expectation  of  the 
interference  parameter  r(u,v)  equals, 

E{r(u,v)}  - 2(p2-l  + p"1)  (4.1) 

which  is  practically  equal  to  the  mathematical  expectation  of  E{r(u,v)} 
for  random  sequences  as  derived  in  Section  2.3.  Furthermore,  the  results 
in  Appendix  D enable  us  to  derive  in  Section  4.3,  the  first  three  central 

A A 

moments  of  the  odd  correlation  functions  Q~  ~(X)  and  8~(2). 

U f V u 

4.3.  Moments  of  the  odd  correlation  functions 

With  the  odd  cross-correlation  function 

®u,vU)  * Cu,v^‘P)  ' OSiSp-1  (4.2) 


one  obtains  with  (D3), 
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Similarly  with  (D4), 

e{$~(X)}  - P_1(p-2^),  l * 0 . (4.4) 

The  variance  of  0-  ~(4)  follows  from  (D5)  and  (D9) 
vartL_,(i)}  = Ele^-W)  - E2{e-_-a)} 

3 ptl  + P 3[(p-2£)2-p] +p  3(p-24)}  (4.5) 

while  with  (D6)  and  (DIO)  one  obtains 

var{0~(4)}  = 4£p  2(p  + l)(p-i),  4 i 0 . (4.6) 

The  third  central  moment  of  0~  ~(4)  is  denoted  as 

u,v 

4—U)  -EC[6.j5(i)  - Et§5-(«})3} 

With  (4.3)  and  (4.5)  this  reduces  to 

~(j t)  = e{q~  ~(i)3  4-  ?(X) 
u,vv  u,v'  ' i 

where 

?(X)  - p'6[(2jl-p)3(2-3p2)  - 3p3  (p2-l)  (2j£-P)3  (4.7) 

and,  after  substitution  of  (D13)  and  (D17), 

El0~}~(i)}  ->  p~2((2-3p)(p-2j&)+6(*)2*  6(P“^)2 +|  i2  (p-i)2  (p-2)  (p-2i)} 

+ 6p"2 (p+1 > { -(f)  (»3  a ) + B3  a ) 3 + ( Pf)  (B3  (P-1 ) + B3  <p-X ) ] 

+ (f  (p-X)  CC^  (i)  +03(1)]  - (^XIC^p-XJ+C^p-X)]}  + 

4-  bp~Z  (p+l)2{** (1)^(1)  + c"(p-X)C3(p-i) 

- 03(4)03 (X)  - B3 (p**X)B^ (p”4)}  . (4.8) 


■v 

„> 

M 


iiVllPrlr  1 : ' 
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For  the  third  central  moment  of  9„(£)  one  obtains  with  (D14), 

(D19),  (4.4)  and  (4.6), 

X~(£)  = e{@~  (4)}  + 6(4) 

where 

6 U)  - p'3(p-22)[4Xa-?)(2  + 3p)  - p2]  (4.9) 

and 

E{9~(4)}  - 8e{c?(A-P)1+12E{C~(4-p)}  + 6EfCu(A-p)}  + 1 

= 2ip‘X{3(2p+l)  -24(3+24)}  + 1 + 48p"1(p+l)B3a)  . (4.10) 

Hence, 

X-(i)  - 8iV3(p+l)[3p-2(p+2)]  + 48p’1(p+l)B3(ja)  . (4.11) 

B^U)  represents  the  number  of  trinomials  (see  Chapter  3)  of 
degree  up  to  but  not  including  4,  which  are  divisible  by  the  primitive 
polynomial  f(x)  which  generates  m-sequence  p,  (or  u). 

represents  the  number  of  trinomials  of  degree  4 and  higher, 
with  the  exponent  of  the  intermediate  term  smaller  than  or  equal  to  4-1, 
which  are  divisible  by  the  primitive  polynomial  f(x)  which  generates 
m-sequcnce  p,  (or  u). 

It  is  shown  in  Appendix  D that 

C3U)  + 3B3 (X)  - C3U)  + 3B3U)  - (2)  . (4.12) 

To  find  33 (4)  --  and  thus  C3U)  --  one  can  use  Lindholm's  equality 

B3U) 

B3  (X)  “IB3U)  - Z d 


(4.13) 
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Here  B ^(i)  equals  the  number  of  trinomials  of  the  form  x^+x^+l, 
with  1 < Cj  < dj  < A- 1,  which  are  divisible  by  f(x). 

Notice  that  (4.8)  and  (4.10)  imply 

- - ^(P-i)  (4.14) 

and 

X-U)  * - X-(p-jt)  (4.15) 

which  reduces  the  number  of  computations  for  each  parameter  with  a factor 

of  1/2.  Furthermore,  for  any  m-sequence  p,  of  length  p = 2n-l  generated 

d,  c , 

by  polynomial  f(x)  of  degree  n,  there  does  not  exist  a trinomial  x J + x J 
+ 1 with  d^  < n,  which  is  divisible  by  f(x).  Hence,  for  l < n 

BjU)  - cn(l)  = 0 , Yu  . (4.16) 

4. A.  Third  moment  evaluation  for  m-sequences  of  length  p = 31  and  p = 63 

In  this  section  the  third  central  moments  of  the  odd  cross- 
correlation and  autocorrelation  functions  as  discussed  in  Section  4.3  are 
evaluated  for  the  m-sequences  of  length  p = 31  and  p = 63.  As  an  example, 
the  values  of  and  C^(^)  for  the  m-sequences  generated  by 

f^(x)  » x^  + x^  + 1 (045),  f ^ (x)  = x^  + x4  -4-  x^  + x^  + 1 (075)  and 
5^2 

f (x)  » x + x + x + x + 1 (067)  respectively,  are  given  in  Table  13. 
Reciprocal  polynomials  give  the  same  values  for  B^(2)  (and  thus  for 

c“(A)). 

To  demonstrate  this  Table  13,  consider  B^(ll)  and  C^(ll)  for 

5 2 

the  m-sequence  u generated  by  f^(x)  = x + x + 1.  The  trinomials  of 

d 4 c , _u  5 2 19 

the  form  x + x J + 1 contributing  to  B^(H)  are  x + x +1  and  x -f 

x*  + 1,  i.e.,  Bj(ll)  “ 2.  ’ 'nee,  “ 2 and  with  (4.12)  one  obtains 
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* • 


Table 

13.  B3CO 

and 

for 

m-sequences  of  length  p = 

31. 

fjtx) 

= X 

5 2 1 
: +x  +1 

f3(x)  = 

5 4 
x +x 

3 2 

■toe  +x  +1 

f5(x)  = 

5 4 
x +x 

2 

+x  +x 

15<x> 

5 3 

f7(x)  = 

5 3 

2 

fn(x)  = 

5 4 

3 

= X 

: +x  +1 

x +x 

+x  +x  +1 

x +x 

+x  +x 

L B*(4)  C*<4) 

l 

B"(i)  C3 W 

l 

B3U)  C3OO 

1 

0 

0 

1 

0 

0 

1 

0 

0 

2 

0 

1 

2 

0 

1 

2 

0 

1 

3 

0 

3 

3 

0 

3 

3 

0 

3 

4 

0 

6 

4 

0 

6 

4 

0 

6 

5 

0 

10 

5 

0 

10 

5 

0 

10 

6 

1 

12 

6 

0 

15 

6 

0 

15 

7 

2 

15 

7 

0 

21 

7 

0 

21 

8 

3 

19 

8 

0 

28 

8 

1 

25 

9 

4 

24 

9 

1 

33 

9 

2 

30 

10 

5 

30 

10 

3 

36 

10 

3 

36 

11 

7 

34 

11 

5 

40 

11 

4 

43 

12 

9 

39 

12 

7 

45 

12 

6 

48 

13 

11 

45 

13 

10 

48 

13 

8 

54 

14 

13 

52 

14 

13 

52 

14 

11 

58 

15 

16 

57 

15 

16 

57 

15 

15 

60 

16 

19 

63 

16 

19 

63 

16 

20 

60 

17 

23 

67 

17 

23 

67 

17 

25 

61 

18 

27 

72 

18 

28 

69 

18 

30 

63 

19 

32 

75 

19 

34 

69 

19 

35 

66 

20 

38 

76 

20 

40 

70 

20 

41 

67 

21 

45 

75 

21 

47 

69 

21 

47 

69 

22 

52 

75 

22 

55 

66 

22 

54 

69 

23 

60 

73 

23 

63 

64 

23 

62 

67 

24 

69 

69 

24 

71 

63 

24 

71 

63 

25 

79 

63 

25 

80 

60 

25 

80 

60 

26 

90 

55 

26 

90 

55 

26 

90 

55 

27 

101 

48 

27 

101 

48 

27 

101 

48 

28 

113 

39 

28 

113 

39 

28 

113 

39 

29 

126 

28 

29 

12t 

28 

29 

126 

28 

30 

140 

15 

30 

140 

15 

30 

140 

15 

31 

155 

0 

31 

155 

0 

31 

155 

0 

52 


s*  . w « 

C^Cll)  ® 34.  The  trinomials  of  the  form  x J + x ^ + 1 and  their 
respective  contributions  to  C^Cll)  are  tabulated  below, 


Trinomial 


5 

, 2 

X 

+ X 

10 

, 4 

X 

+ X 

16 

, 9 

X 

+ X 

18 

X 

+ X 

20 

8 

X 

+ X 

22 

7 

X 

+ X 

27 

, 6 

X 

29 

, 3 

X 

+ X 

+ 1 

+ 1 

+ 1 
+ 1 

+ 1 

+ 1 

+ 1 

+ 1 


Contribution 

3 

6 

2 

10 

3 

4 
4 
2 

Total:  34  • 


The  parameters  X~U)  and  X~  ~(j l)  for  m-sequences  of  length  p = 31 
are  sketched  in  Figure  4 and  Figure  5 respectively.  The  same  parameters 
for  m-sequences  of  length  p ■ 63  are  sketched  in  Figure  6 and  Figure  7. 

4.5.  Discussion  of  X~(X)  and  X~  ~(X) 

■ — u u ,v 

It  is  clear  from  the  results  in  the  previous  section  that  the 
mean  and  the  variance  of  the  odd  correlation  function  values  do  not 
depend  on  the  particular  m-sequences  selected.  The  third  central  moments 
X~(X)  and  X-  ~(X),  however,  clearly  do  depend  on  the  polynomials  f(x)  as 
specified  by  the  expressions  in  Section  4.3,  While  X~(X)  and  X_  -(X),  Vl, 
are  equal  to  zero  for  random  binary  sequences  (see  Section  2.3),  for  the 
m-sequences  they  vary  widely  when  X takes  on  values  in  its  range.  Recall 
that  a third  central  moment  provides  a measure  of  skewness  being  positive, 
zero  or  negative  as  the  distribution  has  a long  positive  tail,  is 


Figure  4 


onm,.* 


Third  central  momenta  o£  the  odd  autocorrelation  function  for 
tn- sequences ; p ■ 31. 


045  (051);  075  (057);  067  (073). 


om  1 1,  j. 


Figure  5 


Hiird  central  moments  of  the  odd  cross -correlation  function 
for  m-sequence3 ; p ■ 31. 


(045,075);  (045,067);  (067,075);  (045,051) 


**-  M'J* 


Figure  7 


Hiird  centre  1 roonients  of  the  odd  cross -correlation  function 
for  m-sequences;  p ■ 63. 


(103,147);  (103,133); 


(147,133) ; (103,141) 
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symmetrical  or  has  a long  negative  tail.  Hence,  the  various  sets 


and 


Xu  3 1 " 1»2»***»P*1} 

Xu,v  3 ^Xu,v^);  ^ 3 O*1****.?-!} 


for  different  m-sequences  and  m-sequence  pairs  respectively  could  be 
potentially  useful  as  sieves  for  sequences  with  good  correlation  properties. 

X * 

For  example,  let  u'  = T u with  x'  chosen  at  random  from  the  integer  values 

A 

in  [0,  p-lj.  The  set  of  values  (0  t (A)  > l = l,2,...,p-l}  for  m-sequence 
u*  can  be  modelled  as  a realization  of  the  set  of  random  variables 

A 

{9 ~ < Jl ) ; l = 1,2, . . . ,p-l} . Hence  one  could  expect  some  correlation  between 

a ^ 

parameters  such  as  0 (u*)  or  min  8 (u1)  and  X~. 

r max v ' , max v u 

x 

On  the  other  hand,  this  correlation  might  not  be  high  enough 

to  show  up  in  the  above  mentioned  parameters  when  the  actual  m-sequences, 

selected  with  X~  as  sieve,  are  compared.  Furthermore,  if  such  a correlation 

exists,  it  is  not  yet  clear  which  measure  on  the  values  in  X-  and  X~  ~ 

u u,v 

should  be  chosen  to  use  as  the  actual  sieve. 

One  measure  might  be  the  presence  of  a certain  number  of  high 

positive  or  negative  peaks  in  X„  or  X-  ~. 

Another  measure  might  be  the  Euclidean  norm  j|X~  J^\  of  the  vector 

(X-  ~(0),  X~  ~(1),...,X~  ~(p-l)),  i.e., 
u,v  u,v  ’ u,v 


M1  • [ X xiv«» 

ImQ 


1/2 


(4.17) 


The  norm  ||xj|  is  defined  in  a similar  way,  with  t i 0.  An  insert  in  the 
Figures  4 through  7 gives  the  Euclidean  norms  for  the  indicated 


m-sequences. 
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A problem  with  both  measures,  and  presumably  will  all,  Is 
the  dominating  and  equalizing  presence  of  the  high  peaks  in  X~  and 
~ for  i < n,  which  do  not  depend  on  f (x). 

4.6.  Actual  data  for  m-sequences  of  length  p = 31  and  p » 63 


In  this  section  we  present  some  actual  data  obtained  for  the 
m-sequences  used  in  Section  4.4.  Additional  data  can  be  found  in 
Sywyk's  work  (1975). 

A 

Let  Fa  denote  the  number  of  times  a certain  value  of  6 (u1) 

max' 

occurs  when  x'  in  u'  = 1*  u takes  on  the  values  0,l,..,,p-l.  Similarly, 

A 

Fc  denotes  the  number  of  times  a certain  value  of  6 (u',v')  occurs 

max 

X * y 1 

when  x'  and  y1  in  u'  * T u and  v'  = TJ  v respectively,  both  take  on  the 

values  0,1,..., p-1.  Fa  and  Fc  are  tabulated  in  Table  14  for  a number  of 

m-sequences  with  the  polynomials  f(x)  in  octal  notation. 

While  a larger  number  of  positive  or  negative  peaks  in 

X~  vor  X~  ~)  or  a higher  value  of.  iiX~ii  (or  j|x~  ~li)  might  indeed  have 

some  positive  correlation  with  a larger  number  of  high  values  of 

8 (u‘)  (or  8 (u',v')),  the  effect  as  a whole  seems  rather  weak, 

max  max 

Actually,  the  data  shows  that  the  m-sequences  do  behave  rather  much  alike. 
Other  results,  collected  while  obtaining  auto-optimal  Gold  sequences  and 
Kasami  sequences  (see  Chapter  5)  do  not  give  more  conclusive  infor- 
mation. 


Hence,  wc  conclude  that  X-  and  X-  while  interesting  on  their 

own  merits,  are  less  useful  as  sieves  for  the  selection  of  o-sequcnces 

with  good  correlation  parameters  such  as  a low  worst -case  value  of 

e <U*)  or  8 (u'.v*). 

max  max 


Table  14.  Fa  and  Fc  for  some  m-sequences  of  length  p = 31  and  p = 63 


P = 31 


Poly. 

9 . (u')  : 

max'  ' 

7 

9 

11 

13 

045  (051) 

Fa  : 

4 

18 

9 

0 

075  (o:7) 

Fa  : 

9 

13 

8 

1 

067  (073) 

Fa  : 

4 

16 

9 

2 

Poly,  pairs 

Wu’-V,) : 

7 

9 

11 

13 

15 

17 

19 

21 

(045,057) 

Fc  : 

0 

53 

339 

333 

156 

56 

20 

4 

(045,073) 

Fc  : 

1 

60 

321 

342 

166 

53 

14 

4 

(057,073) 

Fc  : 

0 

39 

324 

336 

176 

66 

18 

2 

(045,051) 

Fc  : 

10 

188 

336 

216 

121 

67 

21 

2 

p = 63 


Poly. 

9 (u1)  : 

max 

11 

13 

15 

17 

19 

21 

103  (141) 

Fa  : 

8 

7 

22 

16 

9 

1 

147  (163) 

Fa  : 

9 

17 

24 

11 

2 

0 

133  (155) 

Fa  : 

9 

22 

18 

12 

2 

0 

i 

i 


f 
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CHAPTER  5 

CORRELATION  PARAMETERS  FOR  SUMS  OF  PAIRS  OF  M-SEQUENCES 

In  the  previous  chapters  the  discussion  has  been  limited  to 
m-sequences  because  those  sequences  have  excellent  periodic  auto- 
correlation properties.  Sets  of  m-sequences  with  good  periodic  cross- 
correlation properties,  however,  have  a small  cardinality  (see  Section 
3.5).  In  this  chapter  we  will  expand  our  discussion  to  larger  sets  of 
potentially  good  signature  sequences  due  to  Gold  (1967)  and  Kasami  (1966). 


5.1.  Gold  sequences 


Theorem  2 in  section  3.5  yields  a pair  of  m-sequences  of  common 
period  p = 2n  - 1 for  which  the  pairwise  cross-correlation  is  bounded  by 
2 L (n+2 ) / 2j  + ^ Here  n ^ o m(Xi  ^ otherwise  ft(x)  Is  not  primitive. 

In  fact,  the  greatest  common  divisor  of  p and  t in  Theorem  2 is 


gcd(p,t) 


I- 


n # 0 mod  4 
n = 0 mod  4 


(5.1) 


as  was  pointed  out  by  Sarwate  (1976). 

Above  result  is  contained  in  the  stronger  and  more  general 

Theorem  3 which  yields  preferred  pairs  of  primitive  polynomials  for 

n odd  and  c ® 1 or  n = 2 mod  4 and  0*2  but  not  for  n a 0 mod  4.  More 

ll 

preferred  polynomial  pairs  can  be  found  when  t 8 2 + 1 in  Theorem  3 is 

replaced  by  2C  - 2 +1,  as  reported  by  Golomb  (1968).  Furthermore,  the 

number  of  polynomials  for  which  this  theorem  holds  can  be  doubled  with  a 
proposition  in  Pursley  (1976).  This  proposition  states  that  polynomial 
pairs  (fx(x),  ft(x))  and  (f^(x),  f^(x))  generate  m-sequence  pairs  with 
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similar  cross -correlation  values  whenever 

qt  =»  2^  mod(2n  - 1),  for  some  i,  0 — i ^ n-1  . (5.2) 

Analytical  results  such  as  above,  as  well  as  experimental 

results  such  as  discussed  in  section  3.4  yield  m-sequence  pairs  (u,v) 

such  that  ® (u,v)  * + 1.  Furthermore,  it  enabled  Gold  (1967) 

max 

and  Kasami  (1966)  to  formulate  the  following  important  result. 

Theorem  4:  Lot  f^(x)  and  ffc(x)  be  a preferred  pair  of  polynomials  of 
degree  n whereby  n ^ 0 mod  4.  The  shift  register  with  the  product 
f^(x)  ft(x)  as  its  characteristic  polynomial  will  generate  a set  of  2n  + l 
distinct  sequences  of  period  p = 2n  - 1 Any  pair  of  sequences  in  this 

I / i A \ /a  J 

set  has  a cross-correlation  function  bounded  by  2 ^ ' +1,  and  any 

sequence  in  the  set  has  an  autocorrelation  function  whose  out-of-phase 

values  are  also  bounded  by  + 1. 

These  sequences  and  all  their  cyclic  shifts  are  usually  referred 

to  as  Gold  sequences.  The  set  of  distinct  sequences  in  Theorem  4 

consists  of  characteristic  m-sequences  P and  v generated  by  f^(x)  and 

1c 

ffc(x)  respectively  and  sequences  of  the  type  # » n > I v,  k a 0,1,...,  p-1, 

whereby  the  addition  is  trodulo  2.  Hence,  the  Gold  sequences  can  also  be 

generated  by  two  shift  registers  with  characteristic  polynomials  f^(x) 

and  ft(x)  respectively  and  one  modulo  2 adder  (see  Dixon  (1976)). 

Clearly  the  periodic  cross -correlation  function  between  m- 

sequences  p and  v in  the  set  of  Gold  sequences  is  three-valued 

because  f^(x)  and  f^(x)  are  a preferred  pair.  The  periodic  croso -correlation 

Ic 

of  m-seqoence  P and  sequence  a p + T v can  be  regarded  as  a special  case 


Let  d„  (oj,0  denote  the  Hamming  distance  between  any  two 
sequences  m and  £,  both  of  length  p,  and  let  (tu)  denote  the  Hamming 
weight  of  sequence  w.  As  before  we  assume  the  m-sequences  p.  (or  u)  and 
v (or  v)  to  be  in  their  characteristic  form.  Then 

6,  (i)  * p - 2d  («u,A) 

WjZ  fl 

= p - 2WR  (CU+A) 

= p - 2W„  (u+A+A+A^v)  • 

ii 


q (k-H-r(i))  , otherwise  . 
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The  choice  of  p.  and  u characteristic  in  m and  £ is  not  essential  for 
equations  (5.3)  and  (5.5)  and  would  merely  constitutes  a change  in  shift 
4 because 


9-  ~(4)  = 9 (4  +y  - x mod  p) , Vi 

u,v  u,v'  7 r ’ 


(5.6) 


where  u = TXu  and  v * T^  ? . 

Given  a value  of  -^-k,  the  trinomial  equations  in  (5.4) 

r (Jt)  i 

(see  Section  3.2)  specify  for  each  4 the  binary  trinomials  x v + x + 1 
and  x8(£)  + ^ + 1 which  should  be  divisible  by  f^(x)  and  ffc(x) 

respectively. 

For  4 = 1,2,....,  p-1,  the  factor  4 - r(4)  in  equation  (5.5) 

takes  on  all  the  values  in  its  range  and  4 # r (4).  This  implies,  among 

others,  that  9 (4)  in  (5.5)  will  again  be  a three  valued  cross- 

u,w  ° 

correlation  function  when  u and  v are  preferred  m-sequences. 

The  factor  s(4)-r(4)  in  equation  (5.3),  however,  will  take  on 

specific  values  in  [0,p-l].  Those  values  will  depend  on  m-k  as  well 

as  the  trinomial  structure  of  the  m-sequences  involved.  With  the 

difference  s(4)-r(4)  determined  for  all  4,  9 (k  + s(4)-r(4))  still  is  a 

u ,v 

function  of  k. 

In  Che  special  case  that  m=k,  equations  (5.4)  imply  that 
s(l)  - r (4)  €T]^  if  l for  some  i and  y.  Then  if  one  adds  a value  k € T) 
tos(4)-r(4),  the  resulting  sum  k + s (4)  - r (4) fe » the  Gauss'  product 
of  cyclotomic  cosets  ^ and  ^ as  defined  in  (3.17). 

Above  facts  explain  why  certain  choices  of  k in  the  sequence 


<ii  “ u + I v can  yield  a peak  parameter  9 (m)  < 9 (u.v)  and  certain 

max  max 


(k*a)  in  sequence  pair  (m®»n  + I v,  C u + T^u)  can  yield  a peak  parameter 


9 (w,z)  < 9 (u,v).  Here  6 (u,v)  is  assumed  to  be  obtained  from  the 

sso-a  Tiax  max 
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tables  of  Gold  and  Kopitzke  (1965),  Golomb's  theorem  (1968)  or  Kasami's 
results  (1966).  Of  course,  other  correlation  parameters  will  depend  on 
(k,m)  too. 

From  the  results  such  as  above,  we  can  formulate; 

Theorem  5:  Let  f^(x)  and  ft(x)  generate  the  set  of  distinct  Gold 
sequences  H,  u and  u]  = u + Tv,  k = 0,1,...,  p-1  as  in  Theorem  4.  Any 
m-sequence  in  & has  a three-valued  cross -correlation  function  with  any 
other  sequence  in  A sequence  pair  of  the  type  (<a  = |i  + Tku,  C e U + T^v) 
will  have,  as  a function  of  (k,m),  a three-valued  or  in  some  cases  a two- 
valued cross -correlation  function.  Furthermore  any  sequence  in  & of  the 
type  m = U + T^u  will  have,  as  a function  of  k,  an  autocorrelation 
function  with  out-of-phase  values  {-1,-1  + } or  £n  some  cases 

{-1,-1  + 2L(n+2)/2j}. 

It  should  be  noted  that  equation  (5.3)  is  not  restricted  to  pre- 
ferred pairs  of  polynomials.  Hence,  our  investigation  --  in  the  next 
sections  — not  only  includes  Gold  sequences  but  many  other  sequences  .of  the 
type  <u  = IA  + T u.  Furthermore,  we  will  study  in  Section  5.5  the  special  case 

n/2 

t a 2 1,  n even,  i.e.,  polynomial  f^(x)  is  not  primitive.  The  next 

section  illustrates  our  approach  with  an  example  for  sequence  length  p 53  15. 

5.3.  Example  for  m-sequences  of  length  p - 15 

Let  u and  u be  the  two  m-sequences  of  length  p s 15  generated 
4 4 3 

by  polynomials  f^(x)  ® x + x + 1 and  f^(x)  =>  x + x +1  respectively. 

As  was  shown  in  Section  3.4 
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S’ 

-1 

if 

i€l)0 

j 

-5 

if 

0 (£)  - < 
u,v 

3 

if 

7 

if 

It  Il5 

-1 

if 

£€T17 

Hence,  9 (u,v)  * 7. 

max 

Case  1:  m = k 

4 

With  U>  generated  by  f^(x)  3 x + x + 1 one  immediately  obtains 

4 

Tn  + T^  + ua0asa  (trivial)  first  trinomial  satisfying  (5.4),  Hie 

other  trinomials  of  the  form  specified  by  (5.4)  are  tabulated  in  Table  15 

for  the  polynomials  f^(x)  as  well  as  f^ (x) . Observe  from  this  table 

that  s (£)  -r(jJ)  mod  15  € Tl0  , ^ and  . Let  Then  the  sum 

k+s(£)  -r(£)  mod  15  €^1)^,  TJ  and  T]^.  This  implies,  with  Table  5, 

that  k + s(i)  -r(£)  mod  15  t \ and  therefore  8 (£)  * 9 (k +s(/)-r(i)) 

#9  CH-)  3 9 (u,v)  = 7.  Hence,  9 (w)  = 5 for  k€Tu  while  8 (w)  3 7 

u,v  5 max  max  / max 

for  k € Tl?- 
Case  2 : m £ k 

Let  Gx  denote  a series  of  values, 

Gx  3 (x*2^  mod  p)  , j * 0,1,...,  n-1  (5.8) 

and  let 

G(x,y)  • ((x*2^  mod  p,  y-2^  mod  p))  , j = 0,l,...,n-l 


(5.9) 
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Table  15.  Trinomial  structure  of  m-sequences  of  length  p = 15. 


i : 1 

(023)  r(4) : 4 

(031)  a(l):  12 

s(i)-rU):  8 
7+s (l)-r(l):  0 

l : 1 

(023)  r(i): 

(031)  3(1): 

s(4)-r(4) : 


Case  1 : m = k 
2 3 4 5 6 7 

8 14  1 10  13  9 

9 4 3 10  8 13 

1 5 2 0 10  4 

8 12  9 7 2 11 

Case  2:  m = k + 

2 3 4 5 6 7 

8 14  1 10  13  9 

4 3 10  8 13  6 

511  4 9 13  0 12 


8 9 10  11  12  13  14  0 

2 7 5 12  11  6 3 == 

6 2 5 14  1 7 11  == 

4 10  0 2 5 1 8 == 

11  2 7 9 12  8 0 == 

1 

8 9 10  11  12  13  14  0 

2 7 5 12  11  6 3 == 

2 5 14  1 7 11  ==  12 

0 13  9 411  5 ==  == 


( ==  denotes  the  exceptional  cases  in  equation  (5.3)  ) 


Table  16.  Sequence  pairs  (w,z)  for  which  9max(w»2)  < 


W 2 

u-T1 

iC 

V , z 

= u*Tmv; 

p = 15. 

o-k 

k 

(k,m) 

G1 

02, 

G3, 

04, 

G(2 ,3) » 

G(3,4), 

0(4,5), 

03 

00, 

G2, 

G7, 

0(0,3), 

G(2 ,5 ) , 

G(7 , 10) , 

05 

02 

0(2,7) 
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denote  a series  of  pairs.  This  notation  can  simply  be  extended  to  triples 
G(x,y,z)  etc. 

Consider  m^k  + l.  Table  15  indicates  the  resulting  trinomials 

of  the  form  specified  in  (5.4)  and  the  resulting  values  of  s(X)  -x(l)  mod  15. 

Observe  that  k + s(X)  - r(£)  mod  15  £ whenever  k = 2,  3 or  4.  Again 

let  (k,m)  indicate  the  sequence  pair  ((u  = M.+T^v,  C = W+  Trav),  then 

9 (l)  a 0 (k+s(£)  f 6 (u,v)  = 7 for  the  sequence  pair 

w,z  u,v'  x max 

(k,m)  = (2,3),  (3,4)  and  (4,5).  In  fact  one  finds  that  a series  of  values 

m-k  = G1  implies  a series  of  values  k = G2,  G3,  and  G4  or,  alternatively, 

a series  of  sequence  pairs  (k,m)  = G(2,3),  G(3,4)  and  G(4,5)  for  which 

9 (w,z)  t 7 (in  fact  < 5).  Table  16  shows  the  results  for  m-k  = Gx, 

max 

x = 1 , 3 and  5 . 

Examination  of  Table  16  reveals  four  triples  >^3)  w*lere 

^ii 

k^  indicates  sequence  =»  p.  + T Jv,  for  which  the  pairwise  peak 
magnitude  of  the  periodic  cross -correlation  function  equals  5 instead  of  7. 
Those  triples  are  G(i,5,ll). 

In  the  next  section  results  for  sequences  of  longer  length  will 
be  discussed. 


5.4.  Sums  of  pairs  of  m-sequences  up  to  length  p 3 255 

The  trinomial  structure  of  the  o-sequences  u and  v and  its 
relation  with  the  periodic  autocorrelation  and  cross -correlation  functions 
of  w » u + T v,  k€ [0,  p - 1)  has  been  investigated  for  sequence  lengths 
up  to  p * 255.  The  autocorrelation  properties  are  the  least  attractive 
aspects  of  the  sequences  m if  compared  with  m-sequences . Hence,  those 
aspects  received  most  of  our  attention.  In  many  cases  the  AO/LSE  phase 
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shifts  W of  sequences  w were  determined  and  the  cross -correlation  values 
of  interesting  subsets  have  been  computed. 

5.4.1.  Sequence  length  p = 31 

The  trinomials  for  the  m-sequences  of  length  p = 31  are  specified 
in  Table  17  and  are  easily  obtained  from  each  other  with  the  decimation 
property  (3.7)  of  m-sequences. 

Except  in  the  case  that  the  sequences  are  generated  by  reciprocal 

polynomials,  the  m-sequences  p and  v will  be  preferred  pairs,  i.e., 

9max(u,v)  **  l+2^+^^  * 9.  With  the  Gauss'  products  of  cyclotomic 

cosets  of  integers  modulo  31  as  specified  by  Golomb  (1968)  one  obtains 

here  Qmax(w)  ® 2^  + 2 ) / 2 J _ ^ _ -j  £or  exactiy  one  value  of  k namely 

k 6 “Hq  a 0 or  id  = u + u (and  all  cyclic  shifts  of  oj)  . Table  18  gives  the 

results  for  the  periodic  cross -cor relation  parameter  ®max(w>2)* 

asterisk  indicates  that  8 (£)  will  be  a two -valued  function  for  indicated 

w,z 

values  of  (k,m). 

A complete  set  of  periodic  autocorrelation  parameters  C)max(w)» 

La  and  M(w)  for  all  values  of  kSTJ  is  3hown  in  Tcble  19.  Ihe  decimation 
property  (3.7)  and  proposition  (5.2)  imply  here  that  pairs  of  columns  of 
k-values  in  this  Cable  can  be  obtained  from  each  other  by  simple  trans- 
formations of  cosets  T)  . 

x 

A 

As  was  the  case  for  m-3equences,  peak  parameter  ^mfly (tf ) will 
again  depend  on  the  cyclic  shift  of  sequence  w.  Furthermore,  recall 
that  M(v)  + M(v)  * 4S(v),  Vv.  Here  M(v)  is  constant  over  k€Tj^,  some  x. 
Hence,  the  sidelobe  energy  can  be  used  as  a sieve  for  lew  values  of  H(v) 
as  long  as  ktH^,  some  x,  and  Definition  2 (Section  3.3)  can  simply  be 


Table  17.  Trinomial  structure  of  m-sequences  of  length  p = 31. 


2 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15  . 

(045) 

r(2) 

18 

5 

29 

10 

2 27 

22 

20 

-16 

4 

19 

23 

14 

13 

24  . 

(075) 

r(2) 

20 

9 

26 

18 

8 

21 

29 

5 

2 

16 

12 

11 

17 

27 

25  . 

(067) 

r (2) 

19 

7 

11 

14 

29 

22 

2 

28 

15 

27 

3 

13 

12 

4 

9 . 

(057) 

r(2) 

12 

24 

8 

17 

28 

16 

9 

3 

7 

25 

30 

1 

27 

18 

21  . 

(073) 

r(2 ) 

13 

26 

23 

21 

7 

15 

5 

11 

25 

14 

8 

30 

1 

10 

6 . 

(051) 

r(2) 

14 

28 

5 

25 

3 

10 

16 

19 

24 

6 23 

20 

30 

1 

22  . 

ble 

18. 

Good 

8 (w)  or  d (w,z) 

max  max 

; w = u*T^ 

V,  Z s 

u-T^v;  p = 31 

u 

V 

6 (u,v) 

max  * 

e (v)  k 

max 

eytxi) 

9max(w,2)  (k,m) 

045 

075 

s3 

9 

7 

[0] 

7* 

G(  5,10) 

045 

067 

e5 

9 

7 

[0] 

9 

all 

045 

051 

?15 

11 

7 

to] 

7 

G(  1,  2) 

045 

057 

B7 

9 

7 

[0] 

9 

all 

045 

073 

iu 

9 

7 

[0] 

7* 

G(  1,  2) 

Table 

19. 

Periodic 

autocorrelation  functions  for  w 

= u »T  v 

with  k € 7) 

X(W) 

; p = 

31. 

u:  045 

045 

045 

045 

u: 

045 

v:  075 

073 

067 

057 

v: 

051 

(w) 

La 

M(v> 

k 

k 

k 

k 

9 (w) 

max 

La 

M(w) 

k 

7 

10 

510 

[0] 

[0] 

(0) 

[0] 

7 

10 

590 

(0) 

9 

2 

766 

[5) 

m 

(7) 

(15) 

9 

4 

654 

(ID 

9 

4 

830 

(ID 

(15) 

(11) 

(7) 

9 

4 

878 

(7) 

9 

6 

798 

(15) 

(?) 

(D 

(11) 

9 

8 

1134 

(5) 

9 

6 

1086 

(7) 

[5] 

(5) 

(3) 

11 

2 

958 

(3) 

9 

8 

1054 

(1) 

(7) 

(3) 

(1) 

11 

2 

1054 

(15) 

9 

10 

1310 

(3) 

(11) 

(15) 

(5) 

1 1 

2 

1150 

(D 
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extended  to  sums  of  pairs  of  m-sequences  as  well. 

For  the  m-sequences  u and  v generated  by  f^(x)  s 045  and  f^(x) 

= 075  respectively,  Table  44  in  Appendix  E gives  for  each  value  of  k the 

cyclic  shift  y for  which  W * T^(u*T^v)  is  an  AO/LSE  sequence.  As  an 

alternative  specification  of  W,  the  table  also  gives  the  AO/LSE  (octal) 

loading  of  the  shift  register  with  characteristic  polynomial  f^(x)f^(x). 

Tables  45  and  46  give  similar  results  for  the  polynomial  pairs  (f^(x), 

f^(x))  and  (fj(x),  f^(x))  respectively.  Notice  that  in  each  table,  for  each 

value  of  k,  the  AO/LSE  cyclic  shift  y yields  9 „(w)  ^ 8__  (w)I 

max  max 

A good  choice  for  a subset  of  the  AO/LSE  sequences  W in  Table  44 

is  the  set  specified  by  (W:k  j supplemented  with  the  AO/LSE 

m-sequences  U and  V.  In  Table  47,  the  peak  parameters  9mav(w),  ^mx(wiz) 

and  0 (w,z)  for  this  subset  are  tabulated.  As  predicted  by  Table  18, 

max 

the  periodic  cross -correlation  parameter  6njax^w,z^  equals  7 instead  of  9 
for  a number  of  sequence  pairs  (v,z).  The  interference  parameter  r(w,z) 
for  this  subset  is  tabulated  in  Table  48. 

Finally  we  point  out  that  the  data  in  the  Tables  44  through  46 
can  be  related  --  via  equation  (5.3)  --  to  the  various  cases  in 
Figure  5 (Section  4.4).  No  specific  correlation,  however,  could  be 
established  between,  say,  norm  *ad  occurrences  of  ^^^(v)- 

5.4.2.  Sequence  length  p * o3 

As  in  Section  5.4.1,  one  obtains  with  the  trinomial  structure 
and  Gauss'  products  of  cyclotoafc  cosets  of  integers  modulo  p 9 63 
(Table  6)  the  results  tabulated  in  Table  20.  Let  J>'  denote  the  set  of 
distinct  Gold  sequences  v and  m « p + Tkv,  with  p and  v generated  by 
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Table  20.  Good  9 (w)  or  9 (w,z)  ; w = u*Tkv,  z = u*Tmv;  p = 63. 

max  max 


u 

V 

9 (u,v) 

max 

9 

max 

(w)  keytx])  9jnax(w,z) 

(k,m) 

J 

. 0(44,45), 

0(45,50), 

G(37,44) , 

103 

147 

17 

15 

[0], [11], [273 

15*  « 

,0(45,54), 

G(22,25), 

G ( 6 1 ,62) , 

1 

[0(59,62), 

0(55,62), 

G(  0,11) 

103 

155 

23 

15 

[73,  [93, [113, 

9 

0(12,21), 

G( 33 , 42 ) 

[153, C313 

G(  0,  1), 

0(17,18), 

G(20,21), 

G ( 32 , 33 ) , 

G(44,45), 

G(47,48), 

0(55,56), 

G(  1,  4), 

G(20,23), 

0(22,25), 

G(26,29), 

G(27,30), 

0(34,37), 

G(  2,  7), 

G(  1 1 , 16) , 

0(13,18), 

G(20,25), 

G(22 ,27) , 

103 

141 

15 

13 

[133, [213 

13 

0(26,31), 

0(31,36), 

0(39,44), 

G(47,52), 

G(58,  0), 

G(  14,21 ) , 

G(  15,22) , 

0(30,37), 

G(39 ,46) , 

G(  6,15), 

G(  7,16), 

0(11,20), 

0(14,23), 

G(20,29), 

G(29,38), 

G(  9,20), 

G( 10,21 ) , 

0(20,31 ) , 

0(36,47), 

0(46,57), 

G(59,  7), 

G(  0,21), 

G(  6,27), 

G(  9,30) 

103 

163 

23 

15 

C33,  [53,  [73, 

9 

0(21,30), 

G ( 42 ,51) 

[113, [273 

G(  0,  1), 

G(27,3 2), 

0(45,54), 

103 

133 

17 

15 

[03,  [13, [273 

15* 

G(  1,  4), 

G(  1,  8), 

0(16,27), 

0(10,17), 

0(17,20), 

G(  5,10) 
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the  preferred  polynomial  pair  (103,147).  The  data  in  Table  20  shows,  foT 

example,  that  a subset  of  sequences  {to;  ^ ^ "^o '^11*^7^  a peak 

parameter  0 (w)  = 2 ^ + 2)/2j  _ ^ Notice  also  that  in  some  sets 

max 

of  sequences  to,  with  p and  v generated  by  a non-preferred  polynomial 

pair,  a considerable  decrease  of  0 (w)  or  8 (w,z)  can  be  achieved  — 

if  compared  with  8^^ (u,v)  — for  certain  choices  of  k or  (k,m).  Table  21 

specifies  certain  combinations  of  (k,m)  — up  to  four  sequences  in  a 

subset  --  for  which  the  periodic  peak  parameters  are  both  equal  or  better 

than  ®max(u»v)  for  a specific  (u,v)  pair. 

A complete  set  of  periodic  autocorrelation  parameters  is  shown 

in  Table  22.  Furthermore,  all  AO/LSE  sequences  V/  = T^Cu-T^v)  were  obtained 

for  the  polynomial  pairs  (103,147)  and  (103,133)  with  the  results  reported 

In  Appendix  E,  Tables  49  and  50.  For  the  AO/LSE  phase  shifts  of  sequences 

uj  in  the  subset  i«u:  kCH  ,T{  ,T} } of  the  values  of  4 (w) , § (w.z) 

0 11  27  max  max 

and  as  well  as  r(w»z)  are  shown  in  Tables  51  and  52  respectively. 

Observe  that  9 (v,z)  3 15,  in  fact  the  function  9 (£)  is  two -valued 

for  a large  number  of  pairs  of  sequences  as  was  predicted  in  Table  20. 

Of  course,  other  requirements  may  lead  to  different  choices  of 
subsets.  For  example,  one  might  want  to  achieve  a lower  average  value  -« 
over  a subset  --  of  Che  interference  parameter  r(v,z).  Tables  S3  and  54 

show  the  peak  correlation  parameters  and  r(w,z)  for  the  subset 

(W:k  €7]  »T^,  Tt^;  ^ 18}.  For  this  subset,  the  average  value  of  r(w,y)  is 

197.  lower  than  for  subset  (U:k  *r  • 

k 

Let  be  the  set  of  distinct  sequences  v and  <b  * u + T -j,  with 

u and  v generated  by  polynomial  pair  (103,1*1).  a low  value  of  9 (v) 

sax 

= 13  --  with  a (£)  in  this  case  a seven  valued  function  --  is  achieved  for 

V 
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Table 

21. 

Good  0 

max 

(w)  and  8 (w,z)  for  selected  sec  combinations; 

tQaX 

w = U'T^v;  z = u*Tmv;  p = 63. 

u 

V 

8 (w) 

max' 

9 (v,z) 

max'  * 

Good  (k,m)  combinations 

103 

147 

15 

15  < 

i 

r 0(44,45) , G(45,50),  G(37,44), 

G(45,54),  G(22,25),  G(  0,11), 
G(44,45,50),  G(  0,11,44),  G(  0,11,25), 
0(45,54,50),  G( 1 1 ,25,27,54) 

103 

155 

15 

15 

G(15,18),  G ( 30 , 5 1 ) » G( 5 1 ,62) , G(47,50), 
G(55,62),  G ( 37 ,44) , G(  9,18) 

103 

141 

15 

13 

G(  0,  1,  4),  G(20 , 9,31),  G( 17 , 18 ,61 ) , 
G(20, 1 1 ,44) , 008,13,47),  G(22,25,15), 
G(22,25,27),  G(44,39,54) 

103  163 

15 

15 

G( 34 , 37 ) , G(45 ,48) , G(  6,11),  G(10,17), 
G(37 ,44) , G(45,54),  GC 12,33) 

103 

133 

15 

15 

G{  0,  1),  G(27 ,32) , G(45,54), 

G(  1,  4),  G(  1,  8),  G( 16,27) , 

f l t Ty  | U,  1}  U y y VJ 

C(  1,45,54),  G(  1,  4,54,27) 


Table  22.  Periodic  autocorrelation  functions  for  w = u*T^v 


with 

k€T] 

1 (1 

X 

[x]);  p 

= 63. 

u 

: 103 

103 

u 

: 103 

103 

u: 

103 

V 

: 147 

133 

V 

: 155 

163 

v: 

141 

> (v) 

sax 

La 

H(v) 

k 

k 8 (w) 

max 

La 

M(v) 

k 

k 

9 (w) 

sax 

U 

M(w) 

k 

15 

Hi 

3198 

Co] 

[0] 

15 

2 

2110 

[7] 

[7] 

13 

6 

3182 

[21] 

— 

-- 

— 

Cn] 

[1] 

15 

2 

2366 

[31] 

[11] 

13 

6 

3502 

[13) 

— 

— 

— 

[27] 

[27] 

15 

2110 

[93 

[27] 

15 

2 

3646 

C31] 

17 

2 

5118 

C3l] 

[5] 

15 

u 

•i 

2622 

[15] 

[33 

15 

2 

3742 

[11] 

17 

4 

2110 

[15] 

[15] 

15 

4 

3390 

[Hi 

[5] 

15 

2 

3902 

[7] 

— 

-- 

— 

[13] 

C31] 

23 

*■* 

c. 

3166 

[51 

[1] 

15 

2 

4062 

[15] 

— 

— 

— 

[23] 

[11] 

23 

2 

3678 

[23] 

131] 

15 

«* 

34  38 

(23i 

17 

8 

4150 

C3i 

C3] 

23 

2 

4126 

[1] 

[133 

15 

4 

3982 

[3] 

— 

— 

— 

[5] 

[13] 

-- 

— 

[13] 

[233 

15 

4 

4302 

[53 

— 

— 

— 

[7] 

(?] 

23 

4 

6686 

[33 

[15] 

15 

4 

4782 

[1] 

17 

12 

62C6 

[1] 

[23; 

23 

4 

7038 

[27] 

[9] 

15 

6 

2526 

[0] 

— 

— 

— 

[21] 

[21] 

23 

12 

10110 

[0] 

[0] 

13 

6 

4606 

[93 

17 

16 

8254 

[9] 

[9] 

23 

14 

9374 

[21] 

[21] 

15 

6 

50 22 

[27] 
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the  subset  [iu:k  €^3,^}  of*".  A low  value  of  G^Cw.z)  = 13  is  achieved 
for  many  pairs  in  the  subset  {u>:k  € \v\j)  of  *' • AO/hSE  phase  shifts  of 
the  sequences  in  subset  {m:k  € ^3*^21^  an^  «W  are  reported  in 
Appendix  E,  Table  55.  The  peak  correlation  parameters  and  r(w,z)  for  the 
subset  {W:k  €T|lla,T)27a}  are  shown  in  Tables  56  and  57. 

It  should  be  not  -i  that  for  sequence  length  p * 63,  lower 
periodic  peak  parameters  can  be  achieved  for  the  slightly  smaller  sets  of 
Kasami  sequences  (Section  5.5). 

5.4.3.  Sequence  lengths  p = 127  and  p = 255 

For  sequence  length  p = 127  it  is  not  possible  to  specify  values 

of  k in  <u  * u + Tkv  such  that  9 : 9 o (u,v) , whenever  u and  v are 

max  max 

a preferred  pair  of  m-sequences . An  ordering,  however,  of  the  values 

k€T|y  as  a function  of  the  cardinality  La  of  the  set  [4:|9w(i)|  * 

can  be  found  for  the  preferred  polynomial  pairs.  Table  23  gives  the  results. 

k 

Let  * be  the  set  of  Gold  sequences  u,  v and  m = u + T u with 

W and  v generated  by  the  preferred  polynomial  pair  (211,217).  A good  cnoice 

for  a set  of  sequences  in  * with  low  values  of  La  might  be  the  subset 

{(u:k  6 T|q,T|., ,T|^}  supplemented  of  course  with  u and  v.  Table  58  in  Appendix  E 

specifies  the  AQ/LSE  phase  shifts  of  the  sequences  in  subset  {u):k 

as  well  as  £u>:k  € T|^} . The  latter  merely  for  comparative  reasons.  For 

the  subset  {W:k  € 'T]q,T|  } supplemented  with  U and  V,  the  peak  correlation 

parameters  and  r(w,z)  are  given  in  Appendix  E,  Tables  59  and  60.  In 

addition  Table  61  provides  the  cardinality  Lc  of  the  set  {A:|9w^(i)| 

a 9 (w.z)}  for  the  sequences  in  above  subset.  In  the  case  that  z 3 u or 

max 

z 33  v,  Lc  33  28  for  k €T|Q  and  Lc  33  27  for  k € . 

Table  24  gives  results  for  non-preferred  polynomial  pairs. 
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Table  23. 


Periodic  autocorrelation  functions  for  w = u*T  v 


with  1 

«c  €T|  ([x]) 

X 

and 

(u,v)  a 

preferred  pair; 

0 

max 

(w)  = 17; 

p = 

127. 

u:  211 

211 

211 

211 

211 

211 

211 

211 

211 

v:  217 

277 

325 

301 

235 

253 

203 

357 

247 

La 

M(w) 

k 

k 

k 

k 

k 

k 

k 

k 

k 

14 

10430 

[0] 

[03 

[03 

[03 

[03 

[0] 

[03 

[0] 

[03 

18 

13374 

[73 

[293 

[433 

20 

13502 

[193 

20 

14398 

[433 

[13 

[213 

[133 

[193 

[293 

[313 

[473 

[313 

22 

14526 

[33 

[93 

[233 

[93 

[433 

[73 

[113 

[133 

[433 

[293 

[133 

[53 

[553 

22 

15422 

[133 

[53 

[273 

[93 

[153 

[633 

24 

15550 

[193 

[73 

[313 

[553 

[193 

[53 

[53 

[13 

[133 

[293 

[273 

[153 

[73 

[313 

24 

16446 

[153 

[33 

[73 

[73 

[553 

[473 

[21] 

[273 

[233 

26 

16574 

[93 

[33 

[153 

[93 

[193 

[193 

[93 

[293 

[633 

[113 

[273 

[233 

[53 

26 

17470 

[293 

[273 

[433 

[21] 

26 

13438 

[13 

[133 

[S3 

[293 

28 

14462 

[233 

[633 

[213 

[93 

[13 

riii 

[473 

[233 

28 

17598 

[433 

li  u 

[153 

[633 

28 

18494 

[293 

[33 

[293 

[273 

30 

15486 

[13 

[53 

[133 

[133 

[313 

[21] 

[633 

[73 

[553 

[13 

[233 

[193 

[153 

[13 

[213 

[313 

[113 

[553 

30 

18622 

[113 

[33 

[33 

32 

16510 

[113 

[213 

[553 

[53 

[633 

[33 

[153 

[33 

[13 

[213 

[553 

[73 

[633 

[233 

[233 

[193 

[133 

[33 

C313 

[193 

[433 

[93 

[633 

[473 

[553 

[473 

34 

17534 

[473 

[313 

[153 

[113 

[53 

[13 

[273 

[273 

[153 

[213 

... 

[473 

36 

18558 

[473 

[553 

[433 

[233 

[473 

[73 

38 

19582 

[53 

[633 

[433 

[113 

[113 

mum 

C313 

[273 

76 


Table  24.  Good  ; w = u*Tkv  and  (u,v)  a 

non-preferred  pair;  p = 127. 


u 

V 

0 (u,v)  k€\  (Cxi) 

max  x 

9 

max 

211 

367 

41 

C03,  [33,  (53,  (73,  (93,(113, 
[193, [213, (293, (313, [553, (633 

23 

211 

313 

41 

[473 

17 

211 

345 

41 

(313 

17 

211 

221 

21 

[01 

19 

211 

361 

41 

[633 

17 

211 

271 

41 

[633 

17 

211 

375 

41 

[03,  [i],  [73, [193, [213, (233, 

23 

[273, [293,(313, [433, [473, [633 

Finally,  Table  25  gives  some  results  for  sequences  of  length 
p = 255  N<"  preferred  polynomial  pairs  exist  for  this  length  and  good 

periodic  correlation  properties  for  p = 255  can  better  be  achieved  with 
the  small  sets  of  Kasami  sequences  (Section  5.5). 

5.5.  Kasaml  sequences 

In  the  previous  sections  we  investigated  subsets  of  sequences 

IX 

to  = p + T v with  p and  v both  m-sequences  of  period  p « 2 - 1 generated 

by  primitive  polynomials  of  degree  n. 

Let  n be  even.  In  this  section  we  consider  a special  class  of 

sequences  again  of  the  type  is  = P + T t?:  p is  an  m-sequence  of  period 

p - 2n  - 1 generated  by  primitive  polynomial  f^(x)  and  d is  a sequence  of 

period  p'  = 2 ' - 1 generated  by  irreducible  polynomial  f (x)  of  degree 

s 
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Table  25.  Good  Q (w)  ; w = u-Tkv;  p = 255. 

TTlftX 


U V 

e (u,v) 
max 

(Ex3) 

X 

8 (w) 

max 

435  551 

63 

Cl],  [7], [133, [15], [21], [23 3, [25 3, 
[27], [31], C37], [39], [45], [473, [51], 
[533  , [61], [873  , [953  , [111], [127] 

33 

435  747 

65 

[93, [273, [873 

33 

435  545 

47 

[513, [1193 

31 

435  543 

63 

[13,  [93, [253, [273, [293, [313, [373, 
C393, [433, [61], [63], [873, [953 

31 

435  455 

31 

All 

31 

435  703 

95 

[33, [113, [373, [61], [633, [913 

33 

435  607 

63 

[173,[853[1193 

17 

435  561 

31 

All 

31 

435  765 

31 

All 

31 

435  717 

63 

[33,  [93, [133, [193, [21J, [293, [313, 
C393, [533, C553, [953, [ill], [1273 

31 

435  651 

47 

[173, [513 

31 

435  615 

65 

[393, [633, [873 

33 

435  537 

63 

[53,  t'73,  [93, [153, [193, [233  , [253  , 
[273, [313, [373, [393, [453, [473, [513, 
[533  , [593  , [633  , [913, [m3, [1273 

33 
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n/2  whereby  s = 2n^  + 1.  One  can  think  of  v as  a concatination  of  2n^  + 1 

n/2 

m-sequences  u'  of  length  p'  * 2 -1  . A shift  register  with  f1(x)fg(x) 

as  characteristic  polynomial  will  generate  2n^2  sequences  of  period  p = 2n  - 1 
■#  n/2 

(M-  and  o>  = p,  + Tv,  k = 0,1,..., 2 - 1).  Those  sequences,  and  their 

cyclic  shifts,  are  referred  to  as  Kasami  sequences.  Kasami  (1966)  specified 

the  distribution  of  the  Hamming  weight  of  those  sequences.  From  this 

weight  distribution  it  follows  that  the  periodic  cross -correlation  9 ...CO 

takes  on  the  value  8 (t)  = 2n^  - 1 for  2n  * + 2n^2  1 values  of  r and 

u,v' 

®u  v(T)  = ” ^ n ^ ^ ^or  2n  * - 2n//2  * - 1 values  of  t,  where  rC  [0,  p-1] . 

Of  course  0^  v(t2)  s 6y  v(t^)  whenever  mod  p1 . Hence, 

we  can  determine  9 (r) , r € T\  by  specifying  9 (t)  for  r €T1  ' whereby 

U jV  X tl  j V X 

T|  ’ denotes  the  cyclotomic  coset  of  integers  modulo  p' . 

Again  m = p + Tku  and  C = U + Tmv,  with  P and  v' , and  thus  v,  in 
their  natural  orientation.  Immediately  we  have 


9 (w)  < 9 (w,z)  < 2n/2  + 1 . 

max  max 


(5.10) 


As  in  Section  5.2  it  follows  that 


r 


8„ 

w,  z 


P 

-p/p* 

-1 


eu>v(k+8(£)-r(i)). 


if  l *»  0 mod  p,  m-k  = 0 mod  p* 

if  l * 0 mod  p,  m-k  f 0 mod  p' 

if  l t 0 mod  p,  i-fm-k  = 0 mod  p* 

otherwise 


(5.11) 


where  functions  s(^)  and  r(£)  are  specified  as  in  equation  (5.4)  and 
p/p ' = 2n//2  + 1.  As  a special  case  one  obtains,  with  u / w, 
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e 

u,w 


(i) 


-p/p* 


|eu>v(k+*-r(i)). 


if  l “ 0 mod  p 
otherwise  . 


(5.12) 


Notice  that  z(0)[  s 2n^  + 1 = 9max(u>v)  whenever  m-k  / 0 mod  p' . 

Therefore  8 (w,z)  < 8 (u,v)  is  impossible  for  any  subset  of  pairs  of 

max  max  r 

Kasami  sequences,  and  our  main  attention  will  focus  on  the  number  of 

times  8 (w,z)  occurs  for  subsets  of  Kasami  sequences, 

max  ^ 


5.5.1.  Sequence  length  p = 63 

Let  f^(x)  a x^  + x + 1 (103)  generate  m-sequence  li  of  period 

p = 2^  -1  *»  63.  Then  f (x)  = x^  + x^  + 1 (015)  with  s **  2^  + 1 = 9 will 

s 

3 

generate  sequence  v of  period  p'  °2  - 1 a 7.  With  the  matrix  X of 

MI 

Table  7 one  obtains  8 ( l ) « 7 for  l € T|  * and  T)  * while  8 (l)  = - 9 

u,v  0 3 u,v 

for  1 . 


Case  1:  m-k  = 0 mod  p1 


With  the  trinomial  structure  of  m-sequence  u and  v1  it  is  easy 

to  show  that  function  k + s(2)-r(£)  mod  p'  in  (5.11)  does  take  on  values  in 

T1  ' , for  all  k.  Thus  8 (w)  = 8 (u,v)  = 9.  However  the  cardinality 

max  max 

La  of  the  set  {4:|8  (2)|  * 9 (w)}  is  a function  of  k and  one  finds 

w max 


La  = 


if  k€l|  ' 
if  k € nQ » , T)3'  . 


(5.13) 


Case  2:  m-k  ^ 0 mod  p1 


Let 


G'(x,y)  - ((x*2^  mod  p1 


> y* 


(5.14) 
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denote  a series  of  pairs  (x,y) . With  the  same  method  as  used  in  Section 
5-3  it  is  easy  to  show  that  for  sequence  pairs  (uj  = p.  + T^ v,£  » p.  + t“v) 
the  cardinality  Lc  equals 


27 

Lc  = <18 


20 


if  (k,m)  = G' (0,1),  G* (1,3) , G'(l,5),  G' (1,6) 
if  (k,m)  = G' (1,2)  (5.15) 

if  (k,m)  = G' (0,3),  G* (3,5)  . 


Observe  that  G' (1,2,4)  are  triples  for  which  La  = 22  and  Lc  = 18. 

In  the  special  case  that  m-sequence  u is  correlated  with  m =*  p.  + one 
finds  Lc  - 28  for  k^’  and  while  Lc  = 27  for  k € T)^ ' . 

Table  26  specifies  the  AO/LSE  phase  shifts  of  above  discussed 
Kasami  sequences.  The  peak  correlation  parameters  and  r(w,z)  are  given 
in  Tables  27  and  28  respectively. 


5.5.2.  Sequence  length  p - 255 

The  results  for  this  sequence  length  are  of  particular  interest 
because  no  preferred  pairs  of  polynomials  can  be  selected  thus  no  sets 
of  Gold  sequences  exist. 

3^32 

Let  f^(x)  = x + x + x + x +1  (435)  generate  m-sequence  U* 

8 A 

of  period  p = 2 - 1 = 255.  Then  f (x)  * x4  + x + 1 (023)  with 

s 

4 

3=2  + 1 * 17  will  generate  sequence  v of  period  p'  = 15.  One  obtains 
9u,v(2)  = 15  for  *€V  and  V 8u,v(/)  * " 17  £or  1 6V’  V and 

V* 

As  in  5.5.1.  one  has  0 (w)  = 9 (w,z)  = 17  for  all  u>  and 

LIUi  A IUflX 

(w,G).  For  sequences  m,  one  finds  for  cardinality  La 
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Table  26. 

AO/LSE  sequences 

U and  W : 

. Ty 

(u*Tkv); 

u:103 

, v:015 

; P = 

6 (w) 

max 

La 

M(w)  k 

Loading 

y 

8 (w) 

max 

La 

S(w) 

M(w) 

9 

24 

3422  0 

3733 

60 

13 

4 

1455 

2398 

fl 

1771 

33 

11 

4 

1395 

2222 

9 

22 

3358  V2 

4506 

60 

11 

8 

1331 

1966 

l* 

0263 

8 

11 

6 

1427 

2350 

3 

5261 

6 

13 

2 

1483 

2510 

9 

24 

3422  5 

1043 

55 

11 

4 

1423 

2270 

6 

6544 

33 

11 

2 

1247 

1566 

1 

62 

62  U 

0206 

1 

1i 

2 

427 

1646 

Table  27.  Correlation  values 

for 

■ the 

i AO/LSE 

sequences  U and 

W 

= Ty(u 

mk 

•T  v 

);  u: 

103, 

v : 0 1 5 ; 

P = 

63. 

k 

0 

1 

2 

4 

3 

5 

6 

U 

0 

13 

19 

19 

21 

23 

17 

15 

21 

1 

9 

11 

23 

17 

17 

15 

23 

21 

2 

9 

9 

11 

19 

19 

17 

17 

15 

4 

9 

9 

9 

11 

17 

23 

21 

21 

3 

9 

9 

9 

9 

13 

21 

21 

19 

5 

9 

9 

9 

9 

9 

11 

21 

17 

6 

9 

9 

9 

9 

9 

9 

11 

21 

U 

9 

9 

9 

9 

9 

9 

9 

11 

Table  28.  Interference  parameter  r(w,z)  for  the  AO/LSE  sequences  U 


and  W = 

Ty(u  ’ 

Tkv); 

u:  103, 

v:015 

; P = 

63. 

k 1 

2 

4 

3 

5 

6 

U 

0 9570 

7402 

8106 

7338 

6846 

7838 

8962 

1 

6518 

7982 

7030 

6842 

8074 

9022 

2 

6950 

6614 

6962 

6242 

6854 

4 

7742 

6642 

8770 

8830 

3 

6002 

6986 

7638 

5 

7126 

6674 

6 

8282 

82 


La  = 


k«V-V 

teiyV 


For  the  sequence  pairs  (tD,C)  one  finds  for  cardinality  Lc 


(5.16) 


102  if  (k,m)  = 


J~G' (0,5) , G* (0,7);  G' (7,14) , G1 (7,13) ; 
[G* (5*7) , G' (5,14),  G' (5,13) , G*(5,ll). 


104  if  (k,m) 

Lc  => 


119  if  (k,m) 


[o' (1,2), 

O' (1,4); 

G'  (3,6) , 

G'  (3,12) ; 

jo’ (1,3), 

O' (1,6), 

G1 (1,12) , 

G* (1,9) . 

O' (0,1), 

O' (0,3); 

G'  (1,5) , 

G*  (1,10) ; 

O' (3,5), 

O’ (3,10); 

O’ (1,7), 

O’ (1,14), 

G* (1,13) 

, G' (1,11) 

G' (3,7), 

O’ (3,14), 

G' (3,13) 

, G' (3,11) 

In  the  special  case  that  m-sequence  u is  correlated  with  m one  obtains 
Lc  *»  119  for  k € TJ^ * , 11  ’ and  T|^’  while  Lc  = 120  for  k 6 Tt ^ * and 

Table  29  specifies  the  Ao/LSE  phase  shifts  of  above  discussed 
sequences  while  Tables  30  and  31  report  the  resulting  peak  correlation 
parameters  and  r(w,z).  Observe  from  Table  30  that  9c  ° 17  and  s 51 
for  the  set  of  AO/LSE  Kasami  sequences  of  length  p 0 255  generated  by 
polynomial  pair  (f^(x)  ,fg  (x) ) = (103,023).  In  contrast,  =»  95  and 
9£  = 81  for  the  set  of  AO/LSE  m-sequences  of  length  p ® 255  as  reported 


in  Table  421 
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Table  29.  AO/LSE  3equenoes  U and  W = T^(u*Tkv);  u=435,  v:023;  p = 


0 (w) 

max 

La 

M(w) 

k 

Loading 

y 

§ (w) 

max 

A 

La 

S(v) 

M(w) 

17 

110 

60606 

0 

2357 

98 

31 

2 

27807 

50622 

f 

1 

2326 

239 

31 

4 

27127 

47774 

J 

2 

3551 

198 

29 

2 

28227 

52174 

17 

112 

60734 

s 

4 

7430 

35 

23 

U 

29179 

55982 

u8 

5636 

132 

31 

2 

28603 

53678 

3 

3710 

21 

29 

6 

26399 

44862 

6 

5044 

220 

29 

4 

27903 

50878 

17 

112 

60734 

12 

0261 

47 

31 

2 

27211 

48110 

9 

0617 

162 

31 

2 

29287 

56414 

17 

110 

60606 

5 

4666 

89 

31 

4 

29407 

57022 

10 

0454 

89 

29 

4 

26543 

45566 

7 

3134 

198 

31 

2 

27079 

47710 

14 

2577 

58 

31 

2 

28755 

54414 

17 

110 

60606 

13 

3416 

211 

27 

8 

28387 

52942 

11 

7577 

147 

31 

2 

27867 

50862 

1 

254 

254 

U 

4136 

236 

25 

6 

9199 

36542 

Table  30.  Correlation  values  for  the  A0/L3E  sequences  U and 
W = Ty(u-Tkv);  u:435,  v:023;  P = 255. 


k 

0 

1 

2 

4 

8 

3 

6 

12 

9 

5 

10 

7 

14 

13 

11 

U 

0 

31 

51 

43 

41 

41 

37 

51 

47 

39 

45 

45 

43 

45 

45 

41 

51 

1 

17 

31 

51 

45 

43 

49 

41 

41 

39 

49 

39 

43 

41 

43 

41 

45 

2 

17 

17 

29 

45 

41 

47 

47 

41 

51 

45 

49 

41 

47 

49 

45 

43 

4 

17 

17 

17 

33 

47 

41 

41 

39 

45 

41 

45 

43 

43 

41 

41 

41 

8 

17 

17 

17 

17 

31 

43 

39 

45 

41 

41 

37 

39 

41 

43 

51 

41 

3 

17 

17 

17 

17 

17 

29 

49 

41 

41 

39 

47 

45 

39 

47 

37 

49 

6 

17 

17 

17 

17 

17 

17 

29 

49 

45 

39 

47 

43 

43 

49 

47 

43 

12 

17 

17 

17 

17 

17 

17 

17 

31 

41 

41 

45 

43 

45 

39 

41 

45 

9 

17 

17 

17 

17 

17 

17 

17 

17 

?1 

39 

45 

43 

41 

49 

45 

43 

5 

17 

17 

17 

17 

17 

17 

17 

17 

17 

31 

39 

49 

49 

49 
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CHAPTER  6 
CONCLUSIONS 

The  asymptotic  behavior  of  the  aperiodic  correlation  parameters  of 
random  binary  sequences  has  beta  studied  and  an  accurate  approximation  of  the 
average  signal -to-noise  ratio  at  the  correlation  receiver  output  was  obtained. 
This  result  is  very  useful  for  preliminary  system  design.  Those  m-sequences 
that  are  best  suited  for  synchronization  as  well  as  multiple -access  were 
reported.  The  relationship  between  the  third  central  moments  of  the  odd 
correlation  functions  and  the  trinomial  structure  of  m-sequences  has  been 
studied.  Actual  data  showed,  however,  that  the  third  central  moment  is  not 
a very  effective  indication  of  good  aperiodic  correlation  properties.  A new 
method  based  on  Gauss'  products  of  cyclotomic  cosets  was  presented  and  applied 
to  obtain  new  subsets  of  sequences  with  better  correlation  properties. 

The  importance  of  the  data  presented  in  this  thesis  stems  from  its 
use  in  the  selection  of  sequences  for  SSMA  systems.  This  can  be  illustrated 
by  a particular  example.  Suppose  that  there  are  K^S  users  for  a SSMA  system 
and  expression  (2.29)  indicates  that  for  the  system  S^/N^,  Che  required 
sequence  length  is  p = 255.  One  should  select  the  eight  AO/X25E 
m-sequences  as  given  in  Table  41  (p.  104)  if  the  synchronization  and  anti- 
multipath  requirements  necessitate  this.  If,  however,  the  peak  cross- 
correlation  parameters  are  of  primary  concern  one  should  select  eight 
AO/LSE  Kasaai  sequences  from  Table  29  (p.  85)  with  k € tIq » ^5 • ^7  11- 

In  the  latter  table  up  to  sixteen  users  can  be  accommodated . 
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APPENDIX  A 

INTERFERENCE  PARAMETER  FOR  BARKER  SEQUENCES 

2 

In  this  Appendix  we  show  that  r(u,v)  * 2(p  + p-l)  in  any  set 
of  Barker  sequences  of  odd  period  p.  For  all  known  Barker  sequences 
of  odd  period  p. 


C ( l ) 
u 


l odd 


l even  . 


(A.l) 


Each  such  a sequence  gives  rise  to  three  others  under  the  transformations 
(Turyn  and  Storer  (1961)) 


f (-l)Ju 


j 


i 


(-Dj+V 


(A. 2a) 


(A. 2b) 
(A. 2c) 


Substitution  of  (A. 2a)  or  (A. 2b)  into  (1.7)  gives 


C (*)  • (-1)  C (-i)  , VA 

u,v  u.v'  * 


hence 


Cu,v<‘>Cu.v<i+1>  ■ • M 


which  ioplies 


p-1 

z c (i)c  (i+u  - o . 

2®l-p  U,V  U,V 


(A.  3) 
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Furthermore  we  have  with  (A. 2c) 


P-1 


P-1 


S Cu  v(i)Cu  v(i+1)  " , 2 Cu(i)Cu(i+1) 
X-l-p  U»V  U,V  i-l-p  u u 


(A. 4) 


Substitution  of  (A.l)  into  (A.4)  implies  that  (A. 3)  also  holds  for  (A. 2c) 
hence,  we  conclude  with  (1.20)  and  (A. 3)  that  for  p odd 


p-l  2 

r(u,v)  - 2 E c„  m . (A. 5) 

i-l-p  u»v 


It  was  pointed  out  in  Furs ley  and  Sacvate  (1976)  that  the  sum  in  (A. 5) 
2 

equals  p + p - 1 for  Barker  sequences. 
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APPENDIX  B 

BOUNDS  ON  THE  APERIODIC  CROSS -CORRELATION  FUNCTION 


The  cumber  of  sequence  pairs  (u,v)  for  which  C (i)  • r,  with 

U,  V 

0 s X s p-1  and  |r|  s p -l,  equals 

“«.P.r)  ■ 2P  • (B-1) 

To  prove  (2.4)  we  follow  the  procedure  of  Moon  and  Moser  (1963).  Let 
f(p,G)  denote  the  number  of  sequence  pairs  such  that 


C (u,v)  « max  |c  (f)|  * G - G(p)  . 

““  0 S * £ p-1  u’v 

Each  sequence  pair  will  be  counted  p times  if  h(£,p,r)  is  summed  over  all 
Jt  and  r such  that  0 Si  S p-l,|r|  £ G.  Therefore, 


p.£(p,c:)  * V . .e  **)(„?:{«,) 


2-0  ir  s G 


or 


Hence, 


p*f(p,G) 


. p*i  % 

s (2G+l)2ip  Z (p-f)  * 

j£“0 

s (2G+i)22p  fp"*  + 2(P-1>^}. 


2’2p  £(p,G)  S (2G+L)[p’3/2  + 2p*l(?-l)^). 

The  right  hand  side  of  (3.2)  decreases  to  sero  for  p - »,  if  G(p) 
whereby  c > 0. 


(8.2) 


> 


100 

Let  f(p,4,H)  denote  the  number  of  sequence  pairs  such  that 
jc^  vCO|  s H » H(p).  Furthermore,  let  f(p,H)  denote  the  number  of  pairs 
such  that 

imx(u,v)  2 H - H(p). 


Clearly, 


f(p,i,H)  - 2P  Zl  . £ 


r i H 


( ?H  ) 


2**™  £ (P’X) 


,p+£ 


,X-p 


where  the  sum  is  over  those  integers  t such  that  |t-%(p-i)]  > %H. 
Hence, 

f(p,A,H)  S 2P*2P  2 exp{ - (4p)_1H2} 
or 

£ f(p„e,H)  * 2p  22p  exp{*  (4p)  1H2}  . 

I'O 

Hence, 

2~2p  f(p,H)  £ 2p  exp{-  (4p)"‘H2}. 


(B.3) 


The  right  hand  side  of  (B.2)  decreases  to  zero  for  p -♦  «®,  if 

H(p)  * (2-f« ' )(plogp)^,e  1 > 0 or  H(p)  * > 0.  This  concludes  the 

proof 

Conclusion:  If  u and  v are  drawn  at  random  from  the  set  of  all  22p  pairs 
then 


Fr 


{ 


log  C (u,v) 
° max 

H logp 


} 


- 1 * * v(p) 


where  lim  y(p) 

P -4  CO 


0. 
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| : • APPENDIX  C 

5 

j CORRELATION  PARAMETERS  OF  aO/LSE  M-SEQUENCES 


Table  32.  A0/LSE  m-sequencas  of  length  p = 31. 


Poly. 

Loading* 

Poly. 

Loading* 

§ (u) 

max 

A 

La 

S(u) 

045 

11001 

051 

01001 

2 

107 

067 

00011 

073 

01101 

7 

2 

123 

075 

11110 

057 

10010 

7 

2 
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Table  33.  Correlation  values  for  AO/LSE  m-sequences;  p = 31. 


0 

0 

0 

0 

0 

0 

4 

6 

7 

5 

7 

5 

5 

7 

5 

1 

3 

7 

045 

7 

11 

15 

15 

15 

19 

067 

9 

7 

15 

15 

19 

15 

075 

9 

9 

7 

19 

15 

15 

051 

11 

9 

9 

7 

11 

15 

073 

9 

11 

9 

9 

7 

15 

057 

9 

9 

11 

9 

9 

7 

Interference  parameter  r(u,v)  for  AO/LSE  m-sequences;  p s 31. 


045 

067 

075 

045 

2382 

1646 

1990 

067 

2318 

1910 

075 

220  6 

Table  34. 


Table  35.  AO/LSE  m-sequences  of  length  p s 63. 


Foly. 

Loading* 

Poly. 

Loading* 

§ (u)  La 

max 

S(u) 

103 

000010 

141 

011111 

11 

2 

427 

133 

110001 

155 

011001 

11 

2 

503 

147 

100011 

163 

110101 

11 

2 

435 

Table  36.  Correlation  values  for  AO/LSE  m-sequences;  p = 63* 


1 

1 

1 

1 

1 

1 

0 

3 

4 

4 

5 

6 

3 

3 

7 

1 

5 

3 

103 

11 

21 

17 

21 

19 

19 

133 

17 

11 

19 

19 

31 

25 

147 

17 

23 

11 

19 

25 

23 

141 

15 

23 

23 

11 

21 

17 

155 

23 

15 

17 

17 

11 

19 

163 

23 

17 

15 

17 

23 

11 

Table  37.  Interferenoo  parameter  r(u,v)  for  AO/LSE  m-sequenoes;  p * 63* 


103 

133 

147 

103 

9574 

7954 

6958 

133 

10006 

8138 

14? 

9414 

p'Tt'e  ns 

Poly. 

Table  38. 

Loading* 

AO/LSE  m-sequences  of  length  p s 127. 

Poly.  Loading* 

S(u) 

211 

0010000 

221 

1001101 

17 

6 

2183 

217 

0000101 

361 

1111111 

15 

12 

2015 

235 

0001100 

271 

1000101 

17 

10 

2283 

247 

0010111 

345 

0110001 

17 

8 

2255 

277 

1110001 

375 

0101010 

19 

4 

2295 

357 

1110010 

367 

0110101 

17 

4 

2563 

323 

1110111 

313 

1000111 

17 

4 

2203 

203 

1101101 

301 

0010010 

17 

4 

2087 

325 

0000101 

253 

1101100 

19 

6 

2483 

Table  39.  Correlation  values  for  AO/LSE  m-sequences;  p - 127 


2 

2 

2 

2 

2 

3 

3 

2 

3 

2 

3 

2 

q 

3 

3 

3 

3 

2 

1 

1 

3 

4 

7 

5 

2 

0 

2 

2 

6 

7 

4 

7 

6 

1 

0 

5 

1 

7 

5 

7 

7 

7 

3 

3 

5 

1 

1 

1 

5 

5 

7 

3 

1 

3 

211 

17 

33 

31 

27 

27 

31 

39 

27 

25 

37 

45 

33 

47 

33 

35 

29 

37 

27 

717 

17 

15 

31 

31 

27 

29 

31 

37 

31 

45 

35 

23 

37 

33 

29 

29 

31 

29 

235 

17 

17 

17 

27 

31 

29 

29 

35 

29 

33 

23 

37 

29 

33 

27 

31 

33 

31 

247 

17 

17 

17 

17 

31 

29 

33 

29 

33 

47 

37 

29 

35 

33 

29 

29 

23 

53 

277 

17 

17 

17 

17 

19 

27 

33 

33 

39 

33 

33 

33 

33 

37 

31 

29 

29 

35 

357 

17 

17 

17 

17 

17 

17 

31 

39 

27 

35 

29 

27 

29 

31 

41 

31 

31 

35 

325 

17 

17 

41 

17 

17 

17 

17 

41 

41 

29 

29 

31 

29 

29 

31 

37 

27 

27 

203 

17 

17 

41 

17 

17 

41 

17 

17 

37 

37 

31 

33 

23 

29 

31 

27 

37 

25 

325 

17 

17 

17 

17 

41 

17 

41 

41 

19 

27 

29 

31 

53 

35 

35 

27 

25 

37 

221 

21 

41 

41 

41 

41 

41 

41 

17 

17 

17 

33 

31 

27 

27 

31 

39 

27 

25 

361 

41 

21 

17 

41 

41 

41 

41 

41 

17 

17 

15 

31 

31 

27 

29 

31 

37 

31 

271 

41 

17 

21 

41 

17 

41 

17 

17 

41 

17 

17 

17 

27 

31 

29 

29 

35 

29 

345 

41 

41 

41 

21 

'41 

41 

17 

17 

41 

17 

17 

17 

17 

31 

29 

33 

29 

33 

375 

41 

41 

17 

41 

21 

17 

41 

41 

17 

17 

17 

17 

17 

19 

27 

33 

33 

39 

367 

r, 

41 

41 

41 

17 

21 

17 

17 

41 

17 

17 

17 

i7 

17 

17 

31 

39 

27 

313 

41 

41 

17 

17 

41 

17 

21 

41 

17 

17 

17 

41 

17 

17 

17 

17 

41 

41 

301 

17 

41 

17 

17 

41 

17 

41 

21 

-7 

17 

17 

41 

17 

17 

41 

17 

17 

37 

253 

17 

17 

41 

41 

17 

41 

17 

17 

21 

17 

17 

17 

17 

41 

17 

41 

41 

19 

Tabie  40.  Interference  parameter  r(u,v)  for  AO/LSE  nj-ae'juenoes;  p 


21* 


*.17  235  247  277  357  323  203 


325 


211 

217 

235 

247 

277 

357 

323 

203 

325 


41214  33622  3272?  32022  29070  34394  32978  ’i486  30250 

40222  30954  31406  31446  2955"  3^666  34118  32066 

4204*  32570  33714  32974  32246  33106  33310 
40326  3 “'05 4 32002  29546  30486  33634 
41294  31146  31938  33934  34202 

41958  31258  30806  30770 

39870  33538  33598 
39902  32442 
42894 


Table  41.  AO/LSE  m-sequences  of  length  p = 255. 


Poly. 

Loading* 

Poly. 

Loading* 

Wu>  u 

S(tt) 

455 

01101111 

551 

00000111 

27 

2 

9607 

453 

01000010 

651 

00011110 

27 

2 

10927 

435 

10000101 

561 

11111110 

25 

6 

9199 

537 

10110100 

765 

00111100 

25 

8 

8883 

545 

01001100 

515 

11011001 

27 

4 

9555 

543 

10110001 

615 

01100010 

27 

2 

9875 

607 

00100111 

703 

00111101 

27 

4 

9323 

717 

00111110 

747 

10000100 

27 

2 

10211 

Table  42.  Correlation  values  for  AO/LSE  m-sequences;  p = 255. 


4 

4 

4 

5 

5 

5 

6 

7 5 6 5 7 

5 

6 

7 

7 

5 

5 

3 

3 

4 

4 

0 

1 5 5 6 6 

1 

1 

0 

4 

5 

3 

5 

7 

5 

3 

7 

7 1115 

5 

5 

3 

7 

455 

27 

51 

51 

47 

43 

61 

49 

45  61  71  61  55 

43 

77 

61 

63 

453 

65 

27 

53 

53 

81 

51 

51 

55  71  55  51  49 

57 

57 

51 

45 

435 

31 

63 

25 

41 

43 

53 

47 

59  61  51  57  45 

53 

49 

61 

51 

537 

63 

31 

63 

25 

47 

45 

49 

73  55  49  45  53 

47 

41 

41 

61 

545 

31 

95 

47 

65 

27 

61 

45 

57  43  57  53  47 

53 

43 

43 

51 

543 

63 

63 

63 

63 

63 

27 

49 

49  77  57  49  41 

43 

55 

51 

51 

607 

63 

63 

63 

65 

47 

31 

27 

47  61  51  61  41 

43 

51 

53 

53 

717 

47 

65 

63 

95 

31 

47 

31 

27  63  45  51  61 

51 

51 

53 

63 

551 

31 

63 

63 

47 

63 

95 

65 

63  27  51  51  47 

43 

61 

49 

45 

651 

63 

31 

47 

63 

63 

31 

47 

63  65  27  53  53 

81 

51 

51 

55 

561 

63 

47 

31 

31 

63 

65 

95 

65  31  63  25  41 

43 

53 

47 

59 

765 

47 

63 

31 

31 

63 

47 

63 

63  63  31  63  25 

47 

45 

49 

73 

515 

63 

63 

63 

63 

31 

65 

63 

63  31  95  47  65 

27 

61 

45 

57 

615 

95 

31 

65 

47 

65 

31 

63 

63  63  63  63  63 

63 

27 

49 

49 

703 

65 

47 

95 

63 

63 

63 

31 

63  63  63  63  65 

47 

31 

27 

47 

747 

63 

63 

65 

63 

63 

63 

63 

31  47  65  63  95 

31 

47 

31 

27 

Table  43.  Interference  parameter  r(u,v)  for  AO/LSE  m-sequences;  p a 255* 


455 

453 

435 

537 

545 

543 

607 

717 

455 

168870 

126358 

126230 

129578 

125194 

129690 

135586 

132554 

453 

171990 

131678 

134322 

133650 

134914 

127178 

129010 

435 

166870 

131314 

133010 

136034 

122954 

126130 

537 

167422 

129702 

126166 

131454 

131430 

545 

165462 

137750 

120510 

131462 

543 

167638 

129198 

136278 

607 

166438 

133950 

717 

168998 
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APPENDIX  D 

MOMENTS  OF  APERIODIC  CORRELATION  FUNCTIONS  OF  M-SEQUENCES 

In  this  Appendix,  various  moments  of  the  aperiodic  cross- 
correlation and  autocorrelation  fur^tions  are  calculated.  Let  u * TXu 
and  v » T'v,  where  u and  v are  m-sequences  in  their  characteristic  form. 
Assume  that  x and  y are  independent  random  variables  both  uniformly 
distributed  over  the  integers  in  the  range  [0,  p-1].  With  the  aperiodic 
cross -correlation  for  m-sequences  u and  v defined  as  in  (1.7),  the  higher 
moments  of  C~  ~(1)  can  be  calculated  in  a straightforward  manner.  Similarly, 
the  moments  of  C^(A),  X + 0 can  be  obtained.  The  Latter  are  identical  to 
Lindholm's  (1968)  moments  of  M-tuples  when  the  substitution  p - |l|  «M  is 
used. 


D.l.  First  moments  of  the  aperiodic  correlation  functions 


The  conditional  expectation  of  C~  ~(X)  given  1-p  JS  X < 0,  equals 

u,v 

P-1+1 


e{C“  ~(X)/1-P  < 1 < 0 } - e{  E u.  _ V..  } 

u,v ' " - ja0  j-X+x  j+y 


P-1+*  ,2/P-l  \ P-1 

.L  p Vi.  “j -tJ\L  W 


j«0  x-0 


y**0 


- p’2(p+l) 


(D.l) 


which  follows  from  the  property  that 
Similarly  it  follows  that 


p-1 

r 

j-o 


•1  for  any  m-sequence  u. 


l ] 
£ < 
t * 

■>■  e 


eIc^~(A)/o  < i < p-i}  - p~2(P-jfc) 
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(0.2) 


Hence,  from  (D.l)  and  (D.2),  for  all  l 


E{C~  ~(X)3  - p"2(p  - Ul)  . 

U,  V 


(D.3) 


For  the  conditional  expectation  of  CV(X),  given  1-p  < i < 0,  we 

have 

p-l+X  , p-1 

E{CL(l)/l-p  < i < 0}  - E p 2 u u . 

U J«0  x-0  J ***  J+x 

From  the  shift*  and-multiply  property  of  m-sequences  (see  Section  3.2) 
whereby  uj_r+x  ■ uj .44*  uj+x>  ^x>  ^or  some  r an^  2^0,  one  obtains  simply 


p-l+X  1 p-1 

E{c-(X)/1-P  < l < o)  - 2 P 2 U 

j-0  x«0  J +X 


-p_1(P  + i) 


Therefore,  for  all  l + 0 


eCc-(2)3  - p“1(|A|  - p) 


(0.4) 


D.2.  Second  moments  of  the  aperiodic  correlation  functions 


The  conditional  expectation  of  C~  ~(1),  given  1-p  S l < 0 

u>  v 


, p-l+X  ? 

Et(4>5«)/1-|.  S i < 0}  - E{  £ (Uj.^  »J+y)2)  + 


p-2+1  p-l+i,p-l  > ,p-l  . 

2p  2 2 l 2 u.  . u , ) ( 2 v4i  v ) 

J-0  m-J+l  V K-0  *•*«•  VO  J+y  '**' 


Again  from  the  shift -and -multiply  property  of  m- sequences 
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,2  , -2  P”2+*  P’1+A 

Et(^~(A)/l-p  < A < 0}  = p + A + 2p  2 2 1 


j-0  m-j+1 


(P  + Jt)[l  + p'2(p  + A-1>] 


Hence,  for  all-  A, 


E(c?-a)}  - (p-  U!)[l+p"2(p-  |a|  -1)]  . 

U,  V 


(D.5) 


For  the  conditional  expectation  of  C~(A),  given  1-p  < A < 0, 


we  have 


2 , p-2+A  p-l+A  p-1 

E{CZU)/1-P  < A < 0}  » p + A + 2p  I 2 2 u 


u.  u 


j-0  m-j+1  x-0  J+*  n** 


With  A i 0 and  some  r + 0 and  s i 0,  on  can  write 


u.  u.  u u , - u u , , Vx  ♦ 

j-A+x  j+x  m-A+x  m+x  r+x  s+x 

Then,  u . u , * u. , , Vx,  for  some  t because  r ^ s.  Therefore 

’ r+x  s+x  t+x 

2 , p-2+A  p-l+A 

E(CT(A)/l-p  < A < oi  - p + A+2p  2 2 (-1) 

u ” j-0  m-j+1 

- (p  + A)  (1  * p~^(p  + A - 1)]  • 

Hence,  for  all  A ^ 0, 

ElC^CA)}  - (p-  Ut)[l-p’l(p-  Ul  -1)3  . 


(D.6) 


■3  T?v*j'n  '-fri 5^. . ^jrjTRp  wis^kip »f*sw  ! 


I! 


U 


108 


D.3.  Expectation  of  the  interference  parameter  r(u»v) 


Consider  the  product  of  C~  ~(X)  and  C~^~(£+l), 


EiCu,v(i)  C5>~(i  + 1)/1-P  < X < - 1} 


p-l+X  p+X  p-1  p-1 

P'  jo  jo  x"0  jo  Vw*  't*  ’’to 


V„  V 


, p-l+X  p-l+X  p+X  /p-1 

.*•  2 1 ■ n « / M 

j-0 


- -2p  ‘ Z 1+P  - Z Z (jA  “j-X+rVX-l+xA  * n’j+y'mfy> 


j«0  m»0  x»0 

jjfei,  jjfm-l 


u_  „ ,.Jfzv^v 


y-0 


Hence , 


E(C~  ~(X)  CL  ~<X+l)/l-p  < X <-!} 

U)  V U j V 


-2 


£p.  + X)  p *(-p  +X  - l) 


For  all  X one  finds 


(p  - X - l)p"*(-p  - X - 2)  0<X<p-l 

Et°5  v(1)  =3  3<i  + l)1  'I  -2 

u*v  u*v  (p  + X)  P ^(-p  + x-l)  1-P<X<-1 


0.7) 


Substitution  of  0*5)  and  0*7)  in  the  expression  of  the  interference 
parameter  r(u,v)  defined  in  (1.20)  gives 


E{r(u,v)}  - 


P’1  ' ,2 


Z (2e(g~  ~(X)}  + E{C-  -(X)  C-  -(X  + l)}] 

X«l-p  u»v  u’v  u,v 


2(p2+|  p-l  + (3pfl)  + 2 PZ  X p’2(X-2p-l) 
3 X-l 


Hence , 


Etr(u,v)}  * 2 (p2  - 1 + p"1)  . 


0.8) 
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Also  of  interest  is  the  product  of  C~  ~(A)  and  C~  ~(A  - p)  whereby 

UjV  » 

0 < A < p-1. 

ECO^-P)  CqjW/O  < p-1}  - 

A-l  P-1 -A  /p-1  \ /P-1  \ 

" P ' En  Uj^+X  U W ' En  Vj+y  'W+V 

j-0  m-0  x-0  J y»0  J 

because,  indeed  j j*  A+m  whenever  0 < j < A-l  and  0 < m < p-l-A. 
Therefore,  with  C~  ~(-p)  * 0, 

E{c5j-(i  - p)  ~ (jfc)/0  < A < p-1}  - p'2A(p-A)  • 

In  a similar  manner,  it  follows  that 

E(C~U-p)  C~(A)/0  < A < p-1]  - p~XA (A-p)  • 


D.4.  Third  moments  of  the  aperiodic  correlation  functions 


Of  most  interest  is  the  third  moment  of  C~  ~(A). 

U,v 

E{d  ~ (A)/l-p  < A < 0 } - 
U>V 


p-l+A  , p-l+A 

- B{  £ (ttf  ,*,) J}  + 3<P  + A)  El  £ U v } 
j-0  j'X  j j-0  J £ J 


p-L+A  p-3+i  p-2+A 

- 3Et  £ u.  }V .}  + 31  El  I £ 
j-0  J j-0  o-j+i 


p-l+A 

* ViViVY*’. 1 


0>.9) 


(D. 10) 


p~2(p  + A)(3(p  + A)  -Z]  + 


"'v*  ^ 


i 

! 

i . 
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.5  p-3+4 
+ 6p  * [ E 
j-0 


p-2+4  p-1+4  p-1 
EES 
m*j+l  n^ttfl  x*G 


u j - 44ocUm-  4+xun  - 4+x 


3 


p-3+4  p-2+4 
l E E 
j*0  o*j+l 


p-1+4  p-1 

E E v,,  v , v , ) . 

™-i  r-o  i*at*°*r 


Let  1-p  < 4 < 0.  Whenever  there  exist  a three-tuple  (j,m,n)  such  that 


uj-4+x  Um-4+x  " Un-4+x’ 


(D.11) 


and 


p-1 


0<j<m<n<  p-1  + 4 (D.12) 

• p.  Whenever  such  a three-tuple  does 


Che  suln  JQ  uj -f*x  Vl+X  Vl+X 

noc  exists,  i.e.,  Viw  - Vl+x>  h-  the  *“ 

P£  u,  , u . u , • -1.  Of  the  frt*')  possible  three-tuple3  (j.m.n) 

n j-i+x  m-A+x  n-A+x  \ 3 / 

x*0 

for  which  (D+12)  holds,  there  are  B^Cp  + i)  three-tuples  rot  which  (D.ll) 
holds  for  m- sequence  u.  Hence, 


E{C?  ~(4)/l-p  < 4 < 0}  - p'2[3(p  + 4)  -2)  (p  + 4)  + 
u,v 

+ 6p*2[-  (P.f) 

+ (p  + l)Bj(p  + 4]  (-  (p.i) 

+ (p  + l)Bj(p  + i)]  . 


For  all  4 one  obtains 

E{C^~(4))  - p’2(3(p-  |4J)2  -2(p-  |4l)  +6(P  * ^ ^ )2  } - 
6p~2 (p  + I)  (P3^1)  (B^Cp  - Ul  ) -^B^Cp  - |4 j)J 

+ 6p*2(p  + l)2B3(p  - |4t)  (p  - 1 4 j ) . (D.13) 


f 


-W-RW.' mxrv&vs  Ai^T*?**K'^-^!TrW"*^’lrlyWj R.r"^ii-.v- JW"H?-;  \-  JH  .WWWJ-I'1-  l-f PJU V;J -*< 
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A similar  procedure  for  the  autocorrelation  function  gives. 


for  X j*  0, 


E{c2(X)}  - p"1(lxi  -p)3+6p'1(p  + 1)b!J(p- Ut) 


u 


(D.IA) 


as  previously  derived  by  Lindholm  (1968), 


D.5.  Expectation  of  a product  of  aperiodic  correlation  functions 


Finally  we  derive  the  expectation  of  the  product  of  CL  -^(X-p) 


and  C~  ~(X). 

U , V 


l-l 


- ,2  f-i;-1  - - 1 

U V . J 


E[ci  ^(i-p)  CL  ~(i)/0  < X < p-l}  - B{  I (V/j)  * Sn  W + 

» » )“U  J J J**U 

t’l  l-l  p-l/P-1 


vi-2  £-1  ,P“1  P-1 

+ 2p’4  Ze£(2  ) E I ( S £ 

j»0  m^j+1  n ml  x*0  Uj-l+xuo-£-txUn-X,ht  j=0  m»j+l  n *£  y«0  vj+yVm+yVn+y  ' 


V 


ut  0 < X < p-1.  Whenever  there  exists  a turee-tuple  (J,m,n)  such  that 


uj-£+x  Um-£+x  Un -i+x  * * 


(D.  15) 


and 


0 ^ j < o£  i-i  ; £ S n £ p-1 


(D.16) 


p-1 

the  sum  Z u 


**°  //N 

-1.  Of  the  ^2J(p-i)  possible  three-tuple  (j,m,n)  for  which  (D.16)  holds, 
there  are  C^CX)  three-tuples  for  which  (D. 15)  holds. 


u . u » p.  In  the  other  cases  this  sum  equals 

j-X+x  ra-X+x  n-i-fx  r 


With  C~  ~(-p)  • 0,  one  obtains 


e[c~  P)  C~  ~ (X)/0<£  < p-1}  - p“2£(p-x)  + 

U I ’ U | * 


+ 2p*2  (-(*)  (P-X)  + <P  ♦ DC3 (X)  ] l-(2)  Cp-X)  ♦ Cp  + DC3  (X) ) 


(D.17) 


*- v , i ja.m_ .^uniaw.!^^.v..4  u,.'..J..<ivi,!JJUJP,,u..-!ij;k'  Ji  sy.jy.i  ^f^ppivip^ppip 
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A similar  procedure  for  the  auto-correlation  gives 


e{c2(£-p)  C~(i)/0  < 1 < p-i}  - ~p”^£(p~£)  + 2p”^ [“(2)  (P”X)  + (p  + IJC^Ci)] • 

(D.18) 


D.6.  Relationship  between  (i)  and  (X) 

For  any  m-sequence  u one  has  the  relation 

0 < £ < p-1.^ 


Cu(i-pj  + cu(i)  - $U(X)  - - 1 , 


Hence,  the  left-hand  side  of  (D,18)  can  be  written  as 


E[c?(i-p)  C~(£)/0  <iS  p-l]  - 2{-d(X-p)/0  < £ < p-1}  -E{<~(2-p)/0<£<p-l}. 


Substituting  (D.  14)  and  (D.6)  results  in 


ElcjjU-p)  C~U)/0  < l < p-l}  - p~l[£(£2-p  + £-l>  -6(p  + l)  B^(£))  . (D.19) 


A cctap arisen  between  (D.18)  and  (D.19'  gives 

SBjU)  + Cj(X)  - (J)  • (D.20) 

Indeed,  this  relationship  can  be  obtained  from  the  conditions  (D.ll),  (D.12) 
and  (D.16)  inmediate ly , as  follows.  Again  0 < 1 S p-l.  Whenever  (D-ll) 
holds  for  Bj(X)  three-tuples,  i.e.,  whereby  0<j<o<n<  £-1,  th’Ut 
will  be  3Bj(l)  three-tuples  (j,a,n)  such  that  (D- 15)  holds  whereby 
0 < j < o < £-1,  0 < o < £-1.  Notice  that  there  are  possible  choices 
for  the  pair  (j,a)  such  that  0 < j < a < £-1  and  thus  one  concludes  with 
condition  (D.16)  that  (D.20)  must  hold. 
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D.7.  Fourth  moments  of  the  aperiodic  correlation  functions 

Equally  straightforward  as  the  derivations  of  the  third  moment 
of  C-  ~(i)  in  Section  D-4,  it  is  not  difficult  to  show  that 


eCc£  -W)}  “ (p-Ul){3(p-Ul)-2+2p'2(3(p-Ul  -4j(p-U|  -1]} 


2*p'2<P'J^)2  - 24?'2(p  + 1)(p‘^)((b“(p-U|)  +bJ(p-|1|)) 


+ 24p'2<p+l)2B"(p-|i|)»"(p-U|)  . 
where  B^(p-Ul)  denotes  the  number  of  four-tuples  (j,m,n,t)  such  that 


u.  u u » u , Vx 

j-x  m-x  n-x  t-x 


0<  j <m<  a < t < p-Uj  -1  . 


Notice  that  the  underlined  part  of  E(c~  ~(f )}  equals  the  fourth  moment  of 
C~  --(f)  for  random  binary  sequences  (see  Chapter  2). 

U t V 

The  fourth  moment  of  C^U),  l t 0,  is  given  by  tindholm  (1968). 
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APPENDIX  E 

CORRELATION  PARAMETERS  OF  SUMS  OF  PAIRS  OF  M-SEQUENCES 


Table  44.  AO/LSE  sequences  W = T^Cu'T^v)?  u:045,  v:075;  p = 31. 


S 

i 

•t  t 


9 (w) 

max' 

La 

M(w) 

k 

Loading 

y 

9 (w) 

max' 

A 

La 

S(w) 

M(w) 

7 

10 

510 

0 

7150 

20 

5 

8 

203 

302 

I 

5 

3446 

4 

9 

2 

323 

526 

10 

1764 

8 

7 

6 

339 

590 

9 

2 

766  < 

20 

4377 

15 

7 

2 

271 

318 

9 

7017 

28 

7 

£ 

291 

398 

J8 

3354 

10 

7 

10 

355 

654 

11 

3007 

30 

5 

8 

291 

334 

22 

7112 

23 

7 

2 

303 

382 

9 

4 

830 

13 

3125 

8 

7 

6 

319 

446 

26 

1602 

16 

5 

12 

311 

414 

21 

0750 

1 

7 

6 

379 

686 

15 

2776 

27 

7 

2 

267 

270 

30 

3440 

9 

7 

2 

291 

366 

9 

6 

798 

29 

6651 

16 

7 

2 

283 

334 

27 

7046 

7 

7 

4 

319 

478 

23 

4503 

3 

7 

6 

355 

622 

7 

4552 

22 

7 

4 

399 

510 

14 

4421 

26 

7 

2 

375 

414 

9 

6 

1086 

28 

5151 

13 

7 

2 

379 

430 

25a 

2221 

10 

5 

10 

359 

350 

b 

1110 

9 

19 

4060 

22 

7 

4 

399 

510 

1 

3631 

19 

9 

2 

363 

398 

2 

6741 

28 

5 

6 

331 

270 

9 

8 

1054 

4 

3043 

14 

5 

8 

335 

286 

8 

1262 

6 

9 

2 

419 

622 

16 

0655 

2 

9 

4 

479 

862 

3 

5272 

10 

7 

10 

475 

590 

6 

1646 

30 

7 

6 

471 

574 

9 

10 

1310 

12 

7072 

17 

7 

4 

439 

446 

24 

7613 

4 

7 

2 

423 

382 

17 

7002 

14 

7 

4 

439 

446 

? * 


f 
* 


«•■■  i 
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Table  45.  AO/LSE  sequences  W s ^(u-T^v);  u:045.  v:067;  p * 31. 


6 (w) 

max 

La 

M(v) 

k 

Loading 

y 

8 (w) 

max 

La 

S (w) 

M(w) 

7 

10 

510 

0 

1723 

13 

7 

2 

235 

430 

• 

7 

6522 

26 

7 

2 

323 

526 

14 

5730 

2 

7 

2 

275 

334 

9 

2 

766  < 

28 

5333 

2 

7 

2 

283 

366 

25 

3427 

17 

7 

2 

307 

462 

19 

1347 

12 

5 

10 

279 

35C 

11 

4427 

6 

7 

2 

291 

334 

22 

3053 

24 

5 

12 

311 

414 

9 

4 

830 

13 

1021 

22 

7 

2 

291 

334 

26 

3044 

26 

7 

2 

323 

462 

21 

4103 

3 

7 

2 

295 

350 

1 

5354 

9 

7 

2 

319 

478 

2 

4641 

30 

5 

10 

287 

350 

9 

6 

798 

4 

4440 

7 

5 

2 

247 

190 

8 

7553 

9 

7 

2 

295 

382 

16 

6774 

15 

7 

4 

311 

446 

5 

2102 

15 

7 

2 

355 

334 

10 

4014 

1 

7 

4 

427 

622 

9 

6 

1086 

20 

1071 

25 

7 

2 

355 

334 

9 

*^0336 

26 

9 

2 

411 

558 

18 

5214 

6 

7 

2 

367 

382 

3 

3202 

7 

7 

2 

327 

254 

6 

7321 

11 

7 

4 

391 

510 

9 

8 

1054 

12 

6634 

23 

9 

2 

415 

606 

24 

3242 

11 

7 

4 

391 

510 

17a 

6367 

0 

9 

2 

415 

606 

b 

4601 

20 

15 

4741 

10 

7 

6 

455 

510 

30 

3235 

19 

7 

2 

419 

366 

9 

10 

1310 

29 

4417 

12 

7 

2 

435 

430 

27 

2224 

26 

7 

2 

447 

478 

23 

54$  1 

5 

7 

2 

451 

494 

1 

* • 
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Table  46.  AO/LSE  sequences  W = T 


1 (w) 

max 

La 

M(w) 

k 

Loading 

7 

10 

510 

0 

6537 

* 

7 

5773 

14 

5624 

9 

2 

766  < 

28 

3551 

25 

3530 

19 

» 

4376 

3 

1401 

6 

5224 

9 

4 

830 

12 

4621 

24 

5424 

17 

1014 

* 

11 

3676 

22 

0227 

9 

6 

798 

13 

6271 

26 

4317 

21 

7615 

1 

7000 

2 

0200 

9 

6 1086 

4 

1201 

8 

5773 

16 

6056 

5 

3602 

10 

6552 

9 

8 1054 

20 

4025 

9 

2024 

18 

6052 

15 

7057 

30 

2312 

9 

O 

IaJ 

o 

29 

3140 

27 

3225 

23 

4334 

(u'T^v);  u:067,  v:075;  p * 31. 


A A A 


y 

e (w) 
max 

La 

S (w) 

M(w) 

14 

7 

4 

263 

542 

5 

7 

6 

351 

638 

1 

5 

4 

243 

206 

15 

5 

8 

275 

334 

13 

7 

2 

279 

350 

14 

5 

6 

251 

238 

12 

5 

10 

295 

350 

21 

7 

6 

335 

510 

8 

7 

2 

299 

366 

7 

9 

2 

375 

670 

5 

9 

2 

315 

430 

27 

7 

2 

275 

■>'  -> 

18 

7 

2 

315 

-o2 

12 

7 

8 

343 

574 

27 

9 

2 

303 

414 

30 

5 

8 

271 

286 

6 

7 

2 

351 

318 

23 

7 

4 

395 

494 

30 

7 

4 

367 

382 

5 

7 

6 

351 

318 

3 

7 

2 

379 

430 

18 

7 

2 

367 

414 

8 

7 

4 

375 

446 

5 

5 

8 

335 

286 

12 

5 

8 

347 

334 

25 

7 

4 

359 

382 

9 

7 

2 

407 

318 

26 

7 

4 

455 

510 

17 

5 

6 

395 

270 

11 

7 

4 

459 

526 

5 

9 

2 

467 

558 
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?;  ‘ • Table  47.  Correlation  values  for  the  AO/LSE  sequences  U,  V and 

| j W s Ty(u-Tkv);  u:045,  v:075;  k€T]Q  and  ^ ; p = 31. 


k 

0 

5 

10 

20 

9 

18 

U 

V 

0 

5 

17 

15 

11 

11 

11 

11 

13 

5 

9 

9 

15 

11 

11 

13 

13 

15 

10 

9 

7 

7 

13 

11 

17 

15 

11 

20 

9 

9 

7 

7 

13 

11 

13 

13 

9 

9 

9 

9 

7 

7 

11 

11 

15 

18 

9 

1 

9 

9 

7 

7 

13 

13 

U 

9 

9 

9 

9 

9 

9 

7 

15 

V 

9 

9 

9 

9 

9 

9 

9 

7 

Table  48.  Interference  parameter  r(w,z)  for  the  AO/LSE  sequences  U,  V and 
W = Ty(u-Tkv);  u:045,  v:075;  k € T)Q  and  ^ ; p = 31. 


k 5 10  20  9 18  U V 


0 

5 

10 

20 

9 

18 

U 


2006  1942  2050  1846  1830  1742  1702 

1854  2578  2470  1918  2190  2142 

1898  2406  1830  1846  1790 

2106  2146  2090  1858 

1582  2206  2054 
1678  1526 
1990 
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* • Table  51.  Correlation  values  for  the  AO/LSE  sequences  W = T^ (u*T  v); 

i ! us  103 » v:  147;  lc  € T1Q t T(11  and  ; p = 63. 

i 

£ * ' 

i 


k 

0 

11 

22 

44 

25 

50 

37 

27 

54 

45 

0 

13 

27 

23 

21 

21 

27 

19 

17 

19 

21 

11 

15 

11 

21 

17 

19 

25 

19 

19 

17 

19 

22 

15 

17 

13 

23 

17 

21 

19 

23 

19 

17 

44 

15 

15 

17 

13 

17 

23 

19 

25 

23 

19 

25 

15 

15 

15 

17 

13 

19 

21 

23 

23 

17 

50 

15 

15 

15 

15 

17 

11 

21 

19 

19 

19 

37 

15 

17 

15 

15 

15 

17 

11 

23 

23 

19 

27 

17 

15 

17 

15 

15 

17 

15 

9 

17 

17 

54 

17 

15 

15 

17 

15 

15 

17 

15 

11 

19 

45 

17 

17 

15 

15 

17 

15 

15 

15 

15 

11 

Table  52.  Interference  parameter  r(w,z)  for  the  AO/LSE  sequences 


W s 

Ty(u*Tkv);  u:103, 

v: 147; 

k€Tl0 

and  712_ 

t 5 P 

= 63. 

k 

11  22 

44 

25 

50 

37 

27 

54 

45 

0 

8602  8122 

9470 

8314 

8270 

8098 

8830 

9194 

9814 

11 

9142 

8794 

9022 

8498 

9278 

8258 

7798 

9442 

22 

10098 

8662 

8762 

8446 

8474 

7798 

8234 

44 

10690 

9014 

8538 

8902 

9370 

9110 

25 

10506 

8814 

9714 

9214 

10850 

50 

9890 

9798 

7914 

8366 

37 

8506 

8830 

7882 

27 

7746 

7974 

54 

8874 

V t 

Table  53.  Correlation  values  for  the  AO/LSE  sequences  W a T^(u*T  v); 
u:  103,  v:  147;  k€T|3,  1121  and  7lgf  k*l8;  p a 63. 


k 

3 

6 

12 

24 

48 

33 

21 

42 

9 

36 

3 

13 

25 

19 

19 

15 

21 

21 

19 

23 

17 

6 

17 

15 

17 

17 

19 

17 

19 

17 

23 

19 

12 

17 

17 

15 

19 

19 

21 

17 

17 

21 

21 

24 

17 

17 

17 

15 

19 

21 

23 

19 

19 

19 

48 

17 

17 

17 

17 

13 

23 

23 

19 

29 

15 

33 

17 

17 

17 

17 

17 

11 

23 

29 

21 

19 

21 

17 

17 

17 

17 

17 

17 

13 

17 

19 

19 

42 

17 

17 

17 

17 

17 

17 

17 

15 

21 

21 

9 

17 

17 

17 

17 

17 

17 

17 

17 

13 

15 

36 

17 

17 

17 

17 

17 

17 

17 

17 

17 

15 

Table  54.  Interference  parameter  r(w,z)  for  the  AO/LSE  sequences 


1 a 

v V 

r (u-t  v) 

; u: 

103,  v: 

147;  l 

c€T13’ 

^21  and  'Hg, 

k *18; 

P = 1 

k 

6 

12 

24 

48 

33 

21 

42 

9 

36 

3 

10886 

6602 

7870 

6358 

8358 

6922 

7186 

11286 

6342 

6 

8162 

6958 

8606 

7022 

7562 

6266 

7462 

6798 

12 

7298 

5738 

7410 

6582 

4534 

7578 

7034 

24 

8526 

7926 

7242 

6050 

10526 

5998 

48 

7598 

7530 

5114 

6894 

7614 

33 

6554 

6266 

7742 

8646 

21 

7214 

6482 

5490 

42 

5474 

6346 

9 

5638 

1 > re*?  j n?f^' .^i i f f.^,\  vj^tsv  ; 
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Table  55.  AO/LSE  sequences  W = T^(u*T^v);  u:103,  v : 141 ; 
k€t12i,  1113,  Tl11>  and  1^;  p ■ 63. 


» (w) 

max 

La 

M(w) 

k 

Loading 

y 

0 (w) 

max 

A 

La 

S(w) 

M(w) 

13 

6 

3182 

21a 

2223 

41 

13 

2 

1339 

2174 

b 

6053 

2 

42a 

6463 

57 

11 

2 

1207 

1646 

b 

1103 

28 

'13a 

3745 

41 

11 

4 

1247 

1486 

b 

0564 

10 

26a 

7514 

29 

13 

2 

1339 

1854 

b 

3520 

9 

52a 

3262 

48 

13 

2 

1375 

1998 

1 b 

7110 

27 

13 

6 

3502  “S 

41a 

5341 

59 

13 

2 

1431 

2222 

| 

b 

7313 

27 

19a 

4571 

37 

13 

2 

1475 

2398 

b 

0545 

8 

38a 

0367 

56 

11 

2 

1283 

1630 

l b 

7172 

33 

11a 

5062 

44 

13 

2 

1551 

2462 

b 

5615 

9 

22a 

3320 

0 

11 

6 

1635 

2798 

b 

0430 

42 

44  a 

4026 

58 

11 

2 

1311 

1502 

b 

1127 

25 

15 

2 

3742 

25a 

4203 

37 

13 

2 

1463 

2110 

b 

1575 

2 

50a 

7246 

56 

11 

2 

1327 

1566 

b 

2103 

21 

37a 

4326 

57 

13 

2 

1467 

2126 

b 

0721 

33 

27  a 

0431 

62 

13 

2 

1775 

2078 

b 

4507 

38 

15  6 5022  54a 

1601 

56 

13 

4 

1863 

2430 

b 

1445 

17 

45a 

3412 

54 

11 

8 

1843 

2350 

b 

6022 

28 

123 


Table  56.  Correlation  values  for  the  AO/LSE  sequences  W 


T^u.^v); 


u: 

103,  v 

: 141 

; k< 

^lla 

u 

27a 

; p 

= 63 

• 

k 

11 

22 

44 

25 

50 

37 

27 

54 

45 

11 

13 

23 

21 

15 

23 

27 

21 

21 

19 

22 

15 

11 

19 

21 

19 

23 

27 

17 

21 

44 

13 

15 

11 

17 

23 

21 

21 

23 

19 

25 

15 

13 

15 

13 

23 

19 

15 

17 

29 

50 

13 

15 

13 

15 

11 

23 

19 

21 

21 

37 

15 

13 

15 

13 

15 

13 

21 

25 

19 

27 

13 

13 

15 

13 

13 

15 

13 

19 

17 

54 

15 

13 

13 

15 

13 

13 

15 

13 

17 

45 

13 

15 

13 

13 

15 

13 

15 

15 

11 

Table  57.  Interference  parameter  r(w,z)  for  the  AO/LSE  sequences 


W = 

(u*T^ 

v);  u: 

103,  ' 

7:141; 

k€%la 

i»  ^273  » P 

= 63. 

k 

22 

44 

25 

50 

37 

27 

54 

45 

11 

8154 

6366 

5638 

6014 

7650 

6110 

9190 

5962 

22 

8730 

7538 

7650 

6590 

7322 

7034 

8022 

44 

8366 

6150 

6738 

9230 

7694 

6666 

25 

7734 

6858 

6766 

9254 

6658 

50 

7970 

6406 

6966 

7946 

37 

7570 

8122 

6934 

27 

8550 

6866 

54 

7946 

eSHWIIS! BSPPB5W 
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v Ic 

Table  58.  AO/LSE  sequences  W = T^u^T  v);  u:211,  v:217; 
k€T)QfT)7  , T13  and  p = 127. 


e (w)  is 

max' 

M(w) 

k 

Loading 

y 

®max' 

(w)  La 

S(w) 

M(w) 

17  14 

10430 

0 

03121 

103 

21 

2 

5447 

11358 

! 

7 

3235? 

78 

19 

2 

5563 

8878 

14 

36341 

29 

19 

2 

5463 

8478 

28 

66225 

11 

19 

2 

5571 

8910 

17  18 

13374  < 

56 

73375 

103 

19 

2 

5651 

9230 

112 

13367 

85 

21 

4 

6535 

12766 

97 

07553 

109 

19 

6 

6307 

11854 

167 

33553 

32 

13 

8 

4815 

5886 

3 

54364 

55 

17 

6 

5839 

8830 

6 

20430 

60 

17 

8 

5863 

8926 

12 

60012 

46 

19 

2 

6291 

10638 

17  22 

14526 

24 

25226 

47 

19 

4 

6235 

10414 

48 

04210 

75 

19 

2 

5695 

8254 

96a 

51060 

55 

21 

2 

6407 

11102 

b 

24430 

54 

65 

11654 

108 

19 

4 

6739 

12430 

5 

60763 

12 

19 

2 

7527 

10526 

10 

12732 

34 

21 

2 

7775 

11518 

20 

10774 

98 

21 

4 

8327 

13726 

17  38  19582  40 

57743 

100 

17 

8 

6987 

8366 

80 

73550 

83 

19 

2 

7747 

11406 

33 

73731 

33 

19 

4 

7823 

11710 

66 

61443 

29 

19 

4 

7279 

9534 
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Table  59.  Correlation  values  for  the  AO/LSE  sequences  U,  V and 
W = Ty(u*Tkv);  u:211,  v:217;  k€Tl0  and  Tly{  p s 127. 


k 

0 

7 

14 

28 

56 

112 

97 

67 

U 

V 

0 

21 

27 

27 

31 

31 

43 

35 

31 

31 

41 

7 

17 

19 

27 

35 

27 

31 

29 

29 

35 

41 

14 

17 

17 

19 

31 

37 

31 

33 

31 

31 

37 

28 

17 

17 

17 

19 

31 

27 

25 

31 

31 

31 

56 

17 

17 

17 

17 

19 

27 

35 

25 

29 

23 

112 

17 

17 

17 

17 

17 

21 

27 

29 

27 

33 

97 

17 

17 

17 

17 

17 

17 

19 

31 

37 

39 

67 

17 

17 

17 

17 

17 

17 

17 

13 

29 

29 

U 

17 

17 

17 

17 

17 

17 

17 

17 

17 

33 

V 

17 

17 

17 

17 

17 

17 

17 

17 

17 

15 

Table  60.  Interference  parameter  r(w,z)  for  the  AO/LSE  sequences  U,  V an< 
W = Ty(u-Tkv):  u:21 1 , v:217;  k€T.n  and  T|„ : p = 127. 


7 14  28 

56 

0 

31746  31078  30090 

29106 

7 

32474  30750 

30030 

14 

35018 

30586 

28 

29022 

56 

112 

97 

67 

U 


112  97  67  U V 

33878  34338  30942  31350  28798 
32082  32942  2 1626  33034  32042 
34670  34978  34998  32158  34166 
31866  32542  31026  32394  33522 
31258  30854  29946  30482  28882 
31450  33206  32982  32918 
34282  33226  35002 
31670  31886 
33622 


Table  61.  Cardinality  Lc  for  (k,m)  of  sequenoe  pairs 
( v 3 u ‘T^V,  z s u-rV);  u:211,  v:217  ; p * 127. 


(k,o) 

Lo 

G( 

0,  7) 

37 

G( 

7,14) 

20 

G{ 

7,28) 

18 

G( 

7,56) 

19 
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